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ABSTRACT 

A  compilation  of  techniques  of  analysis  for  ortho¬ 
gonally  stiffened,  flat,  rectangular  plates,  with 
various  combinations  of  loading  and  support  conditions, 
is  presented.  The  specific  types  of  plates  considered 
are  sandwich,  corrugated,  rib-reinforced,  and  integrally 
stiffened.  Orthotropic  plate  theory  is  used  as  the 
major  form  of  solution  throughout  the  report. 

A  complete  bibliography,  consisting  of  more  than 
three  hundred  references,  is  included. 
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INTRODUCTION 
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Historically,  the  first  example  of  a  stiffened  plate 
was  a  flat  slab,  reinforced  by  attached  unidirectional 
girders.  Later,  stiffening  was  provided  by  transverse, 
as  well  as  longitudinal,  beams.  Such  configurations  are 
still  prevalent  in  the  construction  of  ships  and  bridges. 
Concern  about  saving  weight,  while  retaining  strength, 
in  the  design  of  aircraft  resulted  in  the  development 
of  sandwich  construction,  in  which  two  thin  plates  which 
carry  tensile  and  compressive  loads  surround  a  light¬ 
weight,  but  relatively  thick,  core  designed  to  transmit 
shear.  Corrugated  plates  also  have  been  a  popular  form 
of  construction  for  many  years.  However,  the  most  recent 
developments  have  been  in  the  area  of  integrally  stiffened 
plates  such  as  "waffle"  and  "dimple."  The  principle  upon 
which  all  of  the  latter  are  based  is  the  removal  of 
material  from  the  neutral  axis  of  a  flat  plate  by  rolling, 
pressing  or  punching,  or  upon  the  reduction  of  an  origin¬ 
ally  flat  thick  plate  by  milling  into  a  thin  plate  with 
ribs . 

It  is  interesting  to  note,  from  the  curve  shown 
below,  that  interest  in  orthogonally  stiffened  plates 
has  been  increasing  steadily.  This,  in  itself,  demon¬ 
strates  a  need  for  this  survey  of  work  in  the  field. 
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History  of  Interest  in  Stiffened  Plates 

A  great  deal  of  work  in  the  area  of  structural 
analysis  of  orthogonally  stiffened  plates  has  been  per¬ 
formed  during  the  past  twenty  years.  Many  unique  tech¬ 
niques  of  solution  have  been  presented  and  a  great 
number  of  specific  problems  have  been  solved.  But,  due 
to  such  a  wide  scattering  of  work  in  this  area,  it  has 
been  extremely  difficult  for  the  practicing  engineer  to 
find,  in  a  reasonable  amount  of  time,  solutions  to 
specific  cases.  For  this  reason,  the  writer  felt  that 
there  was  a  definite  need  for  consolidation  of  this 
research  into  a  single  report. 

This  study,  therefore,  represents  a  survey  of  the 
most  common  analytical  techniques  in  the  area  of 
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orthotropic  plates.  Within  the  body  of  the  paper,  methods 
of  solution  for  sandwich,  corrugated,  rib-stiffened,  and 
integrally  stiffened  flat,  rectangular  plates  under  static 
loading  conditions  are  presented.  In  addition,  a  complete 
bibliography  is  provided  for  researchers  who  desire  to 
study  specific  problems  in  more  detail.  Since  the  titles 
of  most  of  the  listed  papers  quite  satisfactorily  explain 
the  contents  of  the  papers,  the  bibliography  should  prove 
to  be  a  useful  tool. 

« 

Basically,  there  are  three  major  techniques  of 
analysis  of  orthogonally  stiffened  flat  plates.  The 
first  method  is  based  upon  the  replacement  of  the  actual 
plate  with  a  grid  system  of  bars.  The  major  flaw  of 
this  method  of  attack  is  the  neglect  of  torsional 

i 

rigidities  resulting  from  interactions  of  the  bars.  A  i 

second  technique  is  based  upon  energy  concepts.  This 
is  rather  straight-forward  but  errors  are  often  intro- 
duced  due  to  a  failure  to  satisfy  all  boundary  conditions. 

The  third  method  of  analysis  is  orthotropic  plate  theory, 
which  provides  continuity  of  the  plate.  Although  the 
latter  method  requires  rather  difficult  solutions  for 

5 

complex  boundary  conditions,  this  technique  is  the  most 
general  and  most  accurate.  For  this  reason,  orthotropic 
plate  theory  is  employed  throughout  most  of  this  report. 
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The  characteristic  equation  is  first  developed  and  it  is 
then  applied  to  various  plate  configurations,  loading 
cases,  and  support  conditions.  Finally,  experimental 
techniques  for  determining  orthogonality  constants  are 
discussed. 

The  writer  wishes  to  express  his  gratitude  to  the 
U.  S.  Navy  Marine  Engineering  Laboratory  for  initially 
arousing  his  interest  in  this  area.  This  report  is,  to 
a  large  extent,  based  upon  work  the  writer  performed  for 
USNMEL  during  the  summers  of  1964  and  1965.  The  writer 
also  wishes  to  thank  Dr.  Michael  C.  Soteriades  of  the 
Catholic  Univerity  of  America  for  his  kind  advice  and 
to  Dr.  Richard  D.  Mathieu  of  the  U.  S.  Naval  Academy 
for  his  encouragement  during  the  preparation  of  this 
paper.  In  addition,  the  writer  wishes  to  thank 
Mrs.  Kathy  Jones  for  typing  the  manuscript. 


PART  1  -  DIFFERENTIAL  EQUATION  OF  DEFORMED  PLATE 

The  derivation  of  the  differential  equation  of  the 
deformed  orthotropic  plate  appears  in  many  of  the  refer¬ 
ences  (1,  6,  7,  37,  40,  41,  58,  97,  etc.)  cited  here. 

In  the  derivation,  it  is  assumed  that:  the  load 
acting  on  the  plate  is  normal  to  the  middle  plane; 
planes  normal  to  middle  surface  of  undeformed  plate 
remain  normal  after  deformation;  loads  are  reacted  in 
shear  and  flexure;  transverse  shear  deformation  can  be 
neglected  (this  constraint  will  later  be  removed  in  a 
refinement  of  the  theory). 

An  element  of  the  plate  (Figure  l)  is  considered. 
The  directions  shown  are  considered  positive  (right-hand 
rule  is  used  for  moment  vectors),  and  the  coordinate 
axes  are  chosen  to  coincide  with  the  principal  axes  of 
orthotropy . 

1.1  Expressions  for  Stress  and  Strain 

In  order  to  evaluate  the  forces  and  moments  of 
Figure  1,  it  is  first  necessary  to  write  the  expres¬ 
sions  for  the  strains  and  the  stresses.  Letting  the 
subscripts  define  directions,  and  using  Hooke1 s  Law: 
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Figure  1 

Element  dx-dy  (positive  directions) 
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where 

€,Y  *  normal  and  shearing  strains,  respectively, 

*  normal  and  shearing  stresses,  respectively, 

*  moduli  of  elasticity  and  rigidity,  respectively, 
M  “  Poissonfs  ratio. 


Solving  these  expressions  simultaneously  gives  the 


stress  equations: 

<x,  =  ~ 

I  ~ 

rr  _  E 


(*«  > 


^ -  (6-y  ■**  ^ 


*T#y  8 


1  “  '^«y^yK 

Xy 


f  (2) 


It  is  generally  convenient  to  write  the  strains  in 
terms  of  simplified  constants  with  matrix  subscripts. 
Then,  the  expressions  for  strains  (Eqs.  ( 1 ) )  become: 
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But,  due  to  symmetry: 
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and  there  remain  four  elastic  constants  of  the  ortho¬ 


tropic  plate. 


1.2  Geometry 

Consider  a  plate  element  deformed  as  shown  in 
Figure  2. 
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Figure  2 

Deformed  Element  dx-dy 


The  slope  of  the  deformed  surface, 
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Tne  curvature  of  the  deformed  surface, 
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where  p  is  the  radius  of  curvature.  The  "twist”  of 
the  surface, 

e«, 
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If  s  is  the  normal  distance  from  the  middle  sur¬ 
face  of  the  plate  to  an  investigated  point,  the  strains 
can  be  written  in  terms  of  the  geometry: 
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1.3  Final  Stress  Expressions 

The  final  stresses  are  evaluated  from  Eq.  (3), 
by  inverting  the  matrix  of  elastic  constants: 
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Substituting  the  results  of  Eqs.  (8)  into  Eq.  (9): 


Thus,  the  stresses, 
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l.U  Bending  and  Twisting  Moments 


Recalling  that  z is  the  normal  distance  from  the 
neutral  surface,  the  stress  couples  can  be  evaluated 
from  Eqs.  ( 11 )  ,  (12),  and  (13): 
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1.5  Rigidity  Constants 


Tie  bending  moment  expressions  can  be  rewritten  in 
a  simplified  form: 
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*  l2(l-^l.yU„-) 

"  flexural  rigidity  in  x-direction, 

D  .-^1 _ 

*  flexural  rigidity  in  y-direction, 

=  torsional  rigidity. 
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It  should  be  noted,  at  this  point,  that  various 
systems  for  indexing  of  these  rigidity  constants  are 
found  in  the  literature.  This  variety  of  nomenclatures 
is  extremely  confusing  due  to  the  fact  that  a  symbol 
may  appear  in  several  technical  papers  with  a  different 
meaning  attached  to  it  by  each  author.  Most  of  the 
confusion  arises  from  the  assignment  of  symbols  to  the 
torsional  rigidity  constant.  .The  symbol  D  has  been 
chosen  here  because  it  does  not  appear  in  the  reviewed 
literature,  and,  it  is  hoped,  will  help  to  avoid 
misunderstandings . 

1.6  Shear  Forces 


With  reference  to  Figure  1,  the  equation  of 
equilibrium  of  forces  in  the  z  -  direction: 

'IQ*  +  +  P  =  o. 
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The  equation  of  equilibrium  of  moments  about  the 
x  -  axis : 
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From  the  equilibrium  of  moments  about  the  y  -  axis: 
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after  substitution  of 
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Substituting  Eqs.  (17),  ( 18 ) ,  and  (19)  into  Eqs.  (24) 
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and  (25),  and  solving  for  the  shears, 
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1.7  Differential  Equation  of  Flexure 

Substitution  of  Eqs.  (24)  and  (25)  makes  possible 
the  writing  of  the  equilibrium  equation  (23)  in  terms 
of  the  moments  and  the  load: 

-  2  +  -alMy  -  -  p  .  (2*}) 

Bx3y  Byv  r 

Finally,  from  the  substitution  of  Eqs.  (17),  (1$)  and 
(19),  the  general  differential  equation  of  flexure  of 
the  orthotropic  plate  is  written: 


^  "a*4  +  [2D*  *M‘i  +  Dy  Gf  ’  f ‘  ^ 

For  simplicity,  the  final  form  of  this  equation  is  given 


as : 
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where 


D.y  =  I*  +  2  {M"1  D«  +  M"  Dr)  •  (v) 

1.8  Consideration  of  Shear  Effects 

It  is  recalled  that  one  of  the  assumptions  of  the 
above  derivation  was  the  negligible  effect  of  trans¬ 
verse  shear.  It  must  be  recognized,  however,  that  such 
an  assumption  is  not  always  permissible.  In  some  plates, 
such  as  sandwich  (see  Part  5),  a  low  stiffness  gives 
great  importance  to  the  contribution  of  shear  to 
deflection. 

For  this  reason,  several  investigators  have  con¬ 
sidered  this  effect.  Reissner  (Refs.  78,  169)  first 
recognized  this  problem,  and  it  has  since  been  applied 
to  orthotropic  plates  by  Crawford  and  Libove  (Ref.  305), 
Girkman  and  Beer  (Ref.  21),  Libove  and  Batdorf  (Ref.  155), 
March  (Ref.  159),  Medwadowski  (Refs.  62,  63),  Suchar 
(Ref.  93),  and  Wang  (Ref.  177),  among  others. 

Using  matrix  notation  consistent  with  that  of 
Eq.  (3),  the  stress-strain  relationship  is  given  by: 
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Due  to  symmetry, 

and  the  elastic  properties  of  the  plate  are  character¬ 
ized  by  nine  elastic  material  constants.  Thus,  the 
effect  of  transverse  shear,  together  with  that  of 
transverse  normal  stress,  is  considered. 

Medwadowski  (Ref.  62)  applies  nonlinear  theory  of 
elasticity  to  the  derivation  of  the  characteristic 
equation  of  the  orthcptropic  plate.  He  formulates  the 
problem  by  using  a  system  of  sixteen  simultaneous 
equations  —  Eq.  (33),  equilibrium  equations,  and 
strain-displacement  relationships.  The  problem  is 
reduced  to  a  two-dimensional  one  by  the  introduction 
of  body-force  resultants,  X"  and  Y",  and  weighted 
displacements,  u,  v,  and  w.  With  the  introduction  of 
the  Airy  stress  function,  ,  Medwadowski  reduces  the 
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number  of  equations  to  the  following  four: 
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Aj R  constant,  depending  on  elastic  properties 
and  thickness  of  plate, 
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P«»Pa*  loads  acting  on  upper  and  lower  faces  of 
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and  the  remaining  symbols  have  their  normal  meanings, 
as  used  throughout  this  paper. 

The  effect  of  an  elastic  foundation  at  face 
z  =  "t/2  is  considered,  and  a  system  of  four  partial 
differential  equations,  governing  the  oscillations  of 
an  orthotropic  plate,  is  written.  The  system  is  then 
linearized  and  reduced  to  a  single  partial  differential 
equation  of  sixth  order  through  the  choice  of 

F(x, y,i)  =  stress  function. 

The  values  of  the  constants  contained  in  this  fundamental 
expression  are  then  given,  and  the  solution  of  the 
equation  is  taken  to  be  of  the  Levy-type. 

When  the  body-force  terms  and  the  effect  of  trans¬ 
verse  normal  stress  are  neglected  and  elastic  foundation 
modulus  is  set  equal  to  zero,  Medwadowski’s  linearized 
equations  which  govern  the  bending  behavior  reduce  to 
the  less  complex  expressions  of  Libove  and  Batdorf 
(Ref.  155). 

Since  the  Libove-Batdorf  paper  uses  a  definition 
of  plate  rigidities  which  is  somewhat  unique,  no  attempt 
is  made  here  to  make  these  conform  to  the  definitions 
used  in  the  previous  derivations.  Therefore,  let 
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analogous  to  the  constants  defined  by  Eqs.  (17)  and  (18). 
The  twisting  rigidity  retains  its  definition  as  given 
by  Eq.  (19) : 
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In  addition,  Poisson’s  ratios  are  defined  in  terms  of 
curvatures : 
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The  refinement  to  the  orthotropic  plate  theory, 
which  was  presented  in  the  first  part  of  this  section, 
is  introduced  in  the  form  of  shear  stiffnesses, 


Q, 
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(40) 
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The  bending  moment  equations  which  originally  contained 
flexure  terms  only  (Eqs.  (17),  (IS),  and  (19))  are  now 
refined  by  the  consideration  of  the  shear  contribution: 
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Substituting  the  above  equations  into  Eqs.  (24)  and 
(25),  and  applying  the  equation  of  vertical  equilibrium: 


three  equations  in  w,  Qx,  and  Q  are  written: 
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These  equations  are  then  solved  for  w,  Q  ,  and  Q 

A 

by  means  of  determinants. 
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1.9  Summary 


In  general,  orthotropic  plate  theory  which  omits 
the  effects  of  transverse  shear  gives  good  results  for 
orthotropic  or  orthogonally  stiffened  flat  plates. 

Due  to  the  complexity  of  the  problem,  it  is  often 
worthwhile  to  ’’trade  off”  a  small  degree  of  accuracy 
for  a  great  amount  of  computational  labor.  However, 
the  analyst  must  take  care  that  these  errors  do  not 
become  significant. 

Summarizing  the  results  of  the  bending  theory  of 
orthotropic  plates,  the  differential  equation  of  flexure 
is  recalled: 


H  +  ?r>  +  p)  Ojw 


where 
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The  stresses,  as  given  by  Eqs.  (11),  (12),  and  (13), 
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where 


The  practical  application  of  orthotropic  plate 
theory  depends  upon  the  solution  of  the  lateral  de¬ 
flection,  w.  This  function  is  found  by  solving  the 
fundamental  equation  (31)  for  a  given  set  of  boundary 
conditions.  Then,  it  is  possible  to  evaluate  stresses, 
strains,  moments,  and  shears.  The  application  of  the 
theory  to  specific  cases  of  bending  is  considered  in 
the  following  section. 
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PART  2  -  METHODS  OF  SOLUTION  -  TRANSVERSE  LOADING 


The  solution  of  the  basic  differential  equation  (31) 
of  the  orthotropic  plate, 


2D 


*y 


D  ^ 

y  2>j' 
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> 


is  accomplished  by  one  of  the  following  methods: 


2 . 1  Navler  Solution 


A  double  Fourier  series  solution  of  the  form 

i  (*.y)  -  ^  s.m  ^  ^  ,  (48) 

m*  i  mi 

was  used  by  a  majority  of  investigators.  Huber  (Refs. 

37  -41)  was  the  first  to  give  his  attention  to  ortho¬ 
gonally  reinforced  plates  and  his  approach  was  based 
on  the  Navier-type  function.  Some  of  the  other  in¬ 
vestigators  to  employ  this  method  of  attack  were  Csonka 
(Ref.  220),  Giencke  (Ref.  239)  >  Heller  (Ref.  309), 
Kaczhowski  (Refs. 4$,  49)  ,  Nowacki  (Ref.  71  ),  Raville 
(Ref.  16S) ,  Robinson  (Ref.  79  ),  Schumann  (Ref.  27S ) , 
Soper  (Ref.  90),  Timoshenko  and  Woinowsky-Krieger 
(Ref.  79),  and  Tolotti  and  Grioli  (Ref.  290). 
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The  Navier  solution  is  simple  and  straight¬ 
forward,  even  for  complex  loading  conditions.  However, 
the  method  becomes  cumbersome  in  application  due  to  its 
relatively  slow  convergence  property.  It  is  most 
undesirable,  from  a  computational  standpoint,  when 
higher  derivatives  of  the  deformation  function,  w, 
are  involved. 

2.2  Levy  Solution 

The  Levy-t/pe  solution  is  based  on  a  single  series 
expression  of  the  form: 

(Hat 

where  Ym  is  dependent  upon  y  and  independent  of  x.  This 
technique,  which  is  generally  more  efficient  than  the 
above  type,  was  employed  by  several  investigators,  among 
them  Ando  (Ref.  209) ,  Cornelius  (Ref.  12),  Hajek  (Ref.  25), 
Lekhnitsky  (Ref.  58),  Schade  (Refs.  273-276),  and 
Timoshenko  and  Woinowsky-Krieger  (Ref.  97). 

2.3  Other  Solutions 

Most  of  the  recent  work  has  been  concerned  with 
improving  upon  the  computational  labor  involved  in 


Y*  SIM 


mjr* 

a 


(4^ 


solving  the  fundamental  equation.  Solutions  have  been 
submitted  in  many  forms : 


2.3.1  Maclaurinfs  Series 


Rajappa  and  Reddy  {Ref.  77)  have  applied  Maclaurin’s 
series  to  the  problem  of  the  simply  supported  rectangular 
plate  and  written  the  deformation  equation: 


w  (*,y) ,  w  (o,o)  |  (*  £  v  y  -§)  '^r  1  <u9 


with  jhe  origin  of  the  coordinate  system  taken  at  the 
center  of  the  plate. 


2.3.2  Affine  Transformation 


Brilla  (Ref.  &)  makes  use  of  the  transformation, 


to  solve  this  problem  which  satisfies  the  condition: 

D  *  Jo  D  • 

*y  '  *  r 

The  condition  of  Eq.  (52)  was  recommended  by  Huber 
(Ref.  40)  for  reinforced  concrete  slabs,  and  is  dis¬ 
cussed  by  Timoshenko  and  Woinowsky-Krieger  (Ref.  97). 
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2.3.3  Perturbation  1  -chrtoue 


Vinson  with  Brul  i  2nd  Hess  (Refs.  34,  d C 1  —  3  '  3  ) 
use  of  the  same  transformation  as  Brills  (Eq.  (clj) 
then  go  on  to  obtain  a  perturb  at.  ion  excans  i  in  .,e 
a  for  plates  with  rigidity  ratios  which  sa’.isiy  h-.. 
condition , 


The  parameter,  a,  measures  the  deviation  frc.v.  this 
dition  and  is  given  by: 


This  expression,  together  with  the  transformation, 
changes  the  basic  differential  equation  (31)  to: 


,.a:<e 


(S  40 


.  Q4Vv 
V  W  -  * 


3**39* 

A  comparison  between  this  method  and  the  standard 


-E-  . 


techniques  is  made  in  Reference  101.  It  is  shown  tr.at 
for  the  case  of  a  uniformly  loaded  rectangular  orth"'- 
tropic  plate  on  simple  supports,  the  perturbation 
solution  gives  results  of  reasonable  accuracy,  while 
requiring  considerably  less  computation  time  than  the 
Navier  and  Levy  techniques. 


2.3.4  Complex  Variables 

Mader  (Ref.  61)  considers  the  special  condition 
specified  by  Eq.  (52).  Here,  the  basic  differential 


equation  ( 3 1 )  is  written: 
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where 


The  solution  has  the  form 
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where 
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Suchar  (Ref.  93)  introduces  the  complex  variables, 


it  “  X  •+  X- 
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and  parameters  suggested  by  Lekhnicky  (Ref.  58), 


k  - 

i  i, 


k  -  1  +  iM* 
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to  rewrite  the  basic  cl '  1 1  crential  equation  of  the 
orthotropic  plate: 

Ik  $-(>- + k  **M\ 


=  o 


(6o) 


in  which 


<*<(x,v)  =  W  (a,*1) 


is  a  real  solution  of  the  characteristic  equation. 

In  addition,  Morkovin  (Ref.  65)  also  employs  a 
complex  variable  approach  to  the  problem. 

2.3.5  Energy  Technique 

A  direct  method  of  calculating  stresses  in  ortho- 

% 

tropic  plates,  without  the  usual  intermediate  step  of 
calculating  deflections,  is  presented  by  Coull  (Ref.  13). 
The  method  of  least  work  is  employed. 

The  total  strain  energy  due  to  bending  is  given  as: 
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where  the  constants, 
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The  plate  considered  is  supported  along,  and  parallel 
to,  the  y-axis,  and  it  is  free  along  the  two  remaining 
edges . 

The  moments  and  forces  in  the  plate  are  expressed 
in  terms  of  a  single  series,  and  two  loading  cases, 
symmetrical  and  antisymmetrical ,  are  utilized.  The 
solution  involves  an  assumption  of  an  n-th  order 
polynomial  for  the  bending  moment  M  .  which  makes  possible 

A 

the  writing  of  a  set  of  n  linear  differential  equations 
with  constant  coefficients. 

2.3.6  Design  Application 

An  interesting  and  useful  set  of  tables  for  the 
design  of  orthotropic  plates  with  various  edge  supports 
and  under  various  loadings  is  given  in  a  dual-language 
(English  and  German)  book  by  Krug  and  Stein  (Ref.  57). 
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The  basic  equation  (31)  is  used, 
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and  the  following  notations  and  transformations  are 


introduced: 
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Using  the  above  notation,  the  differential  equation  of 
the  unloaded  plate  becomes 
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In  order  to  reduce  the  number  of  parameters  necessary 
to  define  a  particular  plate  from  three  to  two  (X  and  <2  ) , 
the  loading  is  referred  to  the  §  ,  ^  —  coordinate  system 
(Fig.  3)  and  the  ratio, 


C66) 


e  » 


is  introduced. 


Figure  3 

Transformed  Plate 

In  the  book,  charts  are  drawn  for  the  values, 

e  .  2.00;  1.25;  l.00;0.80;0.50. 

X  .  0.80;  0.40;  O  . 

The  dimension  L  ,  from  Eqs.  (63)  and  (64)  is  chosen  by  the 
user  of  the  influence  surfaces,  such  that  the  length  of 
the  shorter  side  of  the  transformed  plate  becomes  either 
20,  or  16,  or  10,  corresponding  to  the  given  e  .  Some 
results  of  this  procedure  are  shown  in  the  following 
chapter. 
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PART  3  -  FLEXURAL  BEHAVIOR 


In  this  part,  the  bending  problem,  for  various 
combinations  of  edge  supports  and  loadings ,  is 
considered.  The  chapter  is  subdivided  into  sections 
according  to  the  support  conditions,  with  further  sub¬ 
divisions  made  for  different  conditions  of  loading. 

In  the  accompanying  drawings,  the  following  notation 
is  used  for  edge  supports: 

(c)  =  clamped , 

(f)  =  free, 

[Sj  =  simply  supported. 


3.1  Simple  Supports 


The  basic  differential  equation  (31)  of  the  ortho¬ 
tropic  plate  is  recalled: 

pv  ^9  w  .  op  w  ps  *3  W 

■  tsx4  f  i  ^y4  ^  * 

where  the  loading  term  is  a  function  of  x  and  y, 


p  -  >  (<-7) 

and  it  can  represent  any  type  of  transverse  loading  on 
the  plate.  The  double  Fourier  series  (Navier)  form  of 
this  function  is : 
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where 


m  =  number  of  sinusoidal  half-waves  in  the 
x-direction , 

n  =  number  of  sinusoidal  half-waves  in  the 
y-direction. 

The  total  deflection  of  the  plate  shown  in  Figure  4 
is  then  calculated  as  the  sum  of  the  partial  deflec¬ 
tions  produced  by  the  partial  sinusoidal  loadings  of 

Eq.  (68).  v 

by 


Figure  4 

Simply  Supported  Plate 


■ac 


The  general  expression  for  the  deformation  of  a 
simply  supported  rectangular  orthotropic  plate,  in 
the  Navier  form,  is  given  by: 


An  alternate  form  of  solution  (discussed  in  Part 
2)  is  the  single-series  expression  attributed  to  Levy. 
This  solution  is  based  upon  the  equation, 

C7o) 

M 

in  which  Ym  is  a  function  of  y  alone  and  the  b  -  edges 
(x  =  0,a)  are  simply  supported. 

Both  methods,  in  addition  to  any  of  the  additional 
techniques  described  in  Part  2,  will  be  utilized  here 
whenever  they  are  applicable. 

3.1.1  Uniformly  Distributed  Load  over  Entire  Surface 

The  case  of  a  simply  supported  plate,  acted  upon 
by  a  load,  pQ,  uniformly  distributed  over  its  entire 
surface  is  shown  in  Figure  5  and  considered  here. 


F igure  5 

Simple  Supports  -  Uniform  Load 


The  Navier  solution,  as  given  in  Reference  97,  is 

first  utilized.  To  evaluate  the  coefficient  A  of 

mn 

Eq.  (69),  it  is  found  from  Eq.  (68)  that 
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For  the  uniform  loading, 

-  p6  ,  (7 z) 
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and  the  coefficient  thus  becomes: 
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The  substitution  of  Eq.  (73)  into  Eq.  (69)  yields  the 
deformation  expression  for  a  uniformly  loaded,  simply 
supported  plate: 


Using  the  Lexvy  form  of  solution,  Ando  (Ref.  209) 
writes  the  deformation  equation: 


33 


) 


c »  c-*m  -sHy  cos  TQDiy  , 

ci 


0  a  subscript  indicating  initial  condition. 


N  ♦  CosV^u 


The  same  problem  is  solved  by  Hess  and  Vinson 
(Ref.  34)  using  a  perturbation  technique.  The  basic 
equation  is  written  in  the  form: 

4  A  i^(y  X  ^  N 
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which  has  the  roots , 
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There  are,  thus,  three  cases  of  solution,  which  are 
dependent  upon  the  elastic  properties  of  the  ortho¬ 
tropic  plate: 

CASE  I: 


The  deformation  function  for  this  case  is  given  by: 


(76) 


CASE  IX: 

V  I >J  H 
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The  coordinates  and  elastic  constats  are  taken 
according  to  Figure  5,  and  the  parameter, 

■v  ir» 

x  ■  T  ' 

Using  these  equations  as  models,  design  curves  are 
plotted  using  the  dimensionless  parameters, 

p.t*  nw* 

as  functions  of  the  three  independent  variables, 

?y  ,  ~ 

D,  D. 
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1.2  Uniformly  Distributed  Load  over  Part  of  Surface 


A  simply  supported  plate,  acted  upon  by  a  uniformly 
distributed  load  over  a  portion  of  its  surface  (Fig.  6) 
is  now  considered. 


!' 


Figure  6 

Simple  Supports  -  Partial  Load  (I) 

First,  the  Navier  solution  is  applied.  Recalling 
the  equation  giving  the  loading  coefficient  (Eq.  (63)) 
and  applying  the  proper  conditions: 


Substitution  of  this  into  Eq.  (69)  gives  the  defor¬ 
mation  expression: 


Next,  the  L^vy-type  solution,  as  given  by  Ando 
(Ref.  209)  is  examined.  The  x-coordinate  is  changed 
(Fig.  7)  to  coincide  with  the  center  of  the  plate. 


tr 


Figure  7 

Simple  Supports  -  Partial  Load  (II) 
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Deflection  expressions  are  written  for  three  regions  of 
the  plate: 


CASE  I :  *  y  &  tit : 


CASE  II:  %  ±  y  £  : 


CASE  III:  ^  : 


The  constants  of  the  above  equations  are  evaluated 
from  boundary  values: 
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In  deriving  the  influence  surfaces  of  their  book 
of  tables  (Ref.  57),  Krug  and  Stein  make  use  of  a 
transformation  from  the  real  structure  (Fig.  8a)  to  a 
"new”  structure  (Fig.  8b). 


a)  Real  Structure 


b)  Transformed  Structure 
Figure  $ 

Simple  Supports  -  Partial  Load  (III) 


The  method  of  transformation  is  given  in  Eqs.  (62) 
through  (65),  and  the  bending  moments  in  the  plate  are: 
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Numerical  values  given  by  these  integrals,  for  given 
geometry  and  constants  of  orthotropy,  are  given  in 
the  tables. 

3.1.3  Concentrated  Force 

The  next  problem  to  be  considered  here  is  that  of 
the  simply  supported  plate  under  the  action  of  a  con¬ 
centrated  load,  P. 

The  writer  applied  the  Navier  form  of  solution  to 
this  case  in  Reference  309.  This  is  considered  first 
(Figure  9) . 


The  load  coefficient  for  the  double-series  solution, 
from  Eq.  (68),  is  written: 

A  .  iE  Sl(1  •»!«  "1?  .  (87) 

db  ^  b 

The  deformation  expression  is  derived  from  the  sub¬ 
stitution  of  this  quantity  into  the  general  equation 
(69): 


Ando  (Ref.  209)  solves  the  problem  using  the  single 
series  solution.  Using  a  coordinate  system  chosen  such 
that  the  x-axis  is  located  at  the  center,  rather  than 
^he  edge,  of  the  plate  (Fig.  10),  he  writes  two  expres- 
.  ions  for  the  deformation: 

CASE  I:  vi  Sr  y  £  %  • 
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Figure  10 

Simple  Supports  -  Concentrated  Load  (II) 


The  constants  of  Eqs.  (89)  and  (90)  are  the  following: 


C  .  Hl^L _  f  -  «imh  "’3U'!  -  (co*h  "'ll* 

'  D,n'Tt'(iSx-)  1  a  \  « 


-t 
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1  r>^TT»(A-x^  L  «  v 

•  JAUH  Sim  *tPE£  j 


As  in  the  previous  case,  Krug  and  Stein  (Ref.  57) 
make  use  of  the  transformation  of  coordinates,  as  illus 
trated  in  Figures  11a  and  lib. 
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Using  the  transformations  given  by  Eqs.  (62)  through 
(65),  the  moments  are  given: 


M«  *  fff  /d‘  p^,f)  M • 

'  .(hi 

‘  aiTf  P(8',)  ' 

Obviously,  superposition  is  used  for  the  calcu¬ 
lation  of  deflections  in  the  first  two  solutions,  when 
more  than  one  concentrated  load  is  acting,  and  the 
latter  solution  (Eq.  (91))  can  be  rewritten  in  the  form: 

M,  ^  yg  £  P  (C „  1<)  Mt  ,  *)  ■  M 

3.1.4  Line  Loading 


In  solving  the  problem  of  the  simply  supported 
plate  under  the  action  of  a  linear  load  parallel  to  the 
x-axis  (Figure  12),  Ando  (Ref.  209)  uses  Eqs.  (89)  and 
(90)  to  give  the  deformation  for  the  ranges  ■£  <  y  <  b/2 
and  y  <  »2.  »  respectively.  The  constants  are  evaluated 
as  follows : 
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c, 

c4 


f.  >>* 
2(J?-Vj 
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L2T 
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CoS  w 


[- 


COSH 


"Bf  4.  (s«MH  ^*Z  COTV4  sBVjj  ; 

*22:  *?  -  ^CoftU  hL^t?  "fAMM  *^y|  > 

4.  (*M«  cotm  "By)]f 


and 


The  load,  F  (If  )  is  represented  as  a  single  trigono¬ 
metric  series. 

Should  one  consider  the  special  case  where  the 
load  acts  on  the  x-axis  (?  =  0),  the  constants  simplify 
to: 


c, 

ct 
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L*  y 
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tamh  "iI2y 
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Figure  12 

Simple  Supports  -  Line  Load  (I) 

A  more  general  inclination  of  the  line  load  is 
considered  by  Krug  and  Stein  (Ref.  57).  As  previously 
described  for  other  loadings,  the  coordinates  are  trans¬ 
formed  from  those  of  Figure  13a  to  those  of  Figure  13b, 
according  to  Eqs.  (62)  through  (65). 

The  bending  moment  expressions  are  then  given  for 
the  general  case  of  variable  loading,  J  : 


dltf 


l  (f CM m  , 

,  v  \Z(t*  \  < 


,r  J>4  ('r/  r 

*Jf  ^  <—  *>  <•  jf CM Mr(W4 


a)  Heal  Structure 


b)  Transformed  Structure 


Figure  13 

Simple  Supports  -  Line  Load  (II) 

For  the  special  case  of  a  load  uniformly  distributed 
along  a  line,  o£  course,  *<fr  )  =fJand  it  is  placed  in 
front  of  the  integral  sign. 


The  problem  of  a  simply  supported  plate  under  the 
action  of  a  uniform  moment  distributed  along  the  y  *=  b 
edge  (Fig.  14)  is  considered  by  Vinson  and  Brull  (Ref.  103 ). 
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Figure  14 

Simple  Supported  -  Edge  Moment 
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Using  the  transformations, 

V'Y  4J%t  1 

the  boundary  conditions  are: 


w(x,b)  •  w(o,y)  - 

3w|  .  5w 

"  ^‘Ly)  V 


w(«,y) 


-  O, 

o, 


MC«,S>  •  M. 
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If  the  perturbation  solution  is  limited  to  cases  where 
the  rigidity  ratios  (D  /D  )  and  (D  /D  )  are  nearly 
equal,  the  basic  equation  (31 )  of  the  orthotropic  plate 


is  rewritten: 

34w  + 

'ax4  3x*3y* 

W 

*  D.  ' 

C*7) 

where 

•“‘O-fe) 

u  Z  (\ 

-  Cp.,/p.^ 

V 

Wy/Q,/  ' 

The  solution  is  taken  as  a  series  in  powers  of 


(<*) 


which  makes  Eq.  ( 97)  become: 

''VTK-.Sisk 


C w) 
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where  the  operator, 

_«  t  4.  h 

V  "  3x*  V}' 

Since  the  curvature  vanishes  along  edge  y  =  b,  the 
moment , 

M  (x,t)  -  M  -  -  D,  |k  <*>>  . 

and  by  substituting  Eq.  (98)  into  Eq.  (100), 


&  <*- « 


*  3y* 


(loo) 


*fcsfV+... 

3x 


The  governing  equation  for  the  first  term  of  the 
solution  is 

Q  **$.(*>?)  -  o, 

and  the  boundary  conditions  are: 


0 oZ) 


art 

Sx* 
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all  ed«es 
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y  *  0 
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The  solution  is  taken  in  the  form: 


*  Y^y)  ♦  f<y>M  .  f,o*> 

where  Y.  satisfies  homogeneous  boundary  conditions 
and 

From  previous  conditions,  the  function, 

«o 

■  Y  [(c-,  +  q,*)  eoSH 

Mil 

+  (C3  +  C4 x)  simh  I  Sim  y  , 

where  the  constants, 
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k  •— 


The  second  term  is  taken  in  the  form, 

•  V©„(x)  S|«  ^~y  ' 

Mil 

where  the  functions  e  (x)  are  evaluated  by: 

6m(*)  ■  (c¥  +  cfcx  ■**A0x4  +  Alx')  cosh  a 

+  (^7  *Cg*  +  A1,Xt+A1Xl)  SIMM  ^ 

where  the  constants, 
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and 


(C.X.4-OX. 
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24 

c  c,x.+^«n. 
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c4^ 

24 
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The  evaluation  of  these  terms  is  necessary  before  the 
deformation  equation  (98)  can  be  solved. 

A  more  general  case  of  a  non-uniform  moment  dis¬ 
tributed  along  any  line  parallel  to  the  x-axis  (Fig.  15) 
is  considered  by  Ando  (Ref.  209). 


Simple  Supports  -  Line  Moment 


The  problem  is  solved  by  using  the  solution  of  Part 
3.1.4  for  a  line  loading.  If  a  linear  load,  such  as 
that  of  Figure  12,  acts  on  the  plate,  and  another 
load  of  the  same  intensity,  parallel  and  very  near 
to  it  is  allowed  to  act  simultaneously  in  the  opposite 
direction,  a  line  couple  is  created.  This  is  shown 
in  Figure  16. 


1 


Figure  16 

Simple  Supports  -  Line-Load  Couple 
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Thus,  as 

(v,-v.)  —  o , 

the  loading  of  Figure  16  approaches  that  of  Figure  15, 
and 

MOO  -  F(x)  (Vl.V>)  . 


Thus,  the  deformation  for 

two  regions 

is  again  given 

by  Eqs.  (89)  and  (90),  and  the  constants  for  this 

loading  condition  are: 
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and  the  loading  term, 


E. 


D«'T* 


M  (5) 


% 


tODi  dt  . 
<*  ^ 


The  problem  is  greatly  simplified,  of  course, 
when  the  line  of  the  moment  coincides  with  the 
x-axis.  Then,  the  constants  become: 

r  -  E-.W 

“  Z  (**-$) 

r  E.vy 

*  “  2(>'-'X) 


3.2-2  Opposite  Edges  Clamped  and  the  Other  Two  Simply 
Supported 

Since  the  problem  of  an  orthotropic  rectangular 
plate  having  two  opposite  sides  built  in  and  the  two 
remaining  sides  simply  supported  is  of  some  practical 
importance,  it  deserves  discussion  here.  Several 
loading  conditions  are  considered  and  the  method  of 
Ando  (Ref.  209)  is  utilized. 
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3.2.1  Uniformly  Distributed  Load 


A  plate  with  a  uniformly  distributed  loading 
acting  over  the  entire  length,  b,  and  over  a  portion 
of  the  length,  a,  is  considered  (see  Figure  17). 

The  notation  (§)  is  again  used  to  indicate  simple 
supports,  and  (c)  designates  clamped  edges. 


Figure  17 

-  Uniform  Load 
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The  deformation  equations  (83),  (84),  and  185) 
are  again  used  as  the  basic  expressions.  Due  to  the 
boundary  conditions : 


>  (*0*0 


*r  y  •  - 1  * 


and  because  of  the  simplification  of  the  loading 
(relative  to  that  shown  in  Figure  7),  the  constants 


become : 


cs=c,* 


SlMH 


wTXy 


>CO*H  ' 


_ 

a  d  0 


These  values  are  substituted  into  Eqs.  (83),  (84), 
and  (85)  and  deflections  are  found  in  any  of  the  three 
regions.  For  the  case  of  a  uniform  loading  over  the 
entire  surface,  §1  =  0  and  ^  =  a,  in  the  evaluation 
of  the  load  parameter, 
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ZaAf- 
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3.2,2  Concentrated  Force 

Next,  the  orthotropic  rectangular  plate  of  Figure 
10  is  considered.  Here,  however,  the  edge  conditions 
are : 

x  =  0,  a  - -  simply  supported, 

y  =  -  b/2  - *•  clamped. 

Eqs.  (39)  and  (90)  are  used  to  describe  the  deformed 
condition,  and  the  constants  for  these  boundary 
conditions  become: 
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In  ,the  special  case  of  the  load  P  acting  on  the 


x-axis , 

*2  -  o  > 


and  the  constants  become: 

C  .  — 1 -  -  c.  . 
1  E>  *•  IT  (x-V)  a 


When  the  concentrated  load  is  located  at  the  center 
of  the  plate, 


=>  0 


and  the  constants  are  simplified  further: 

C  -  *  f^L  (\ Y  -  AY  \  m  C 


-  (-0*'  p«m 
'*  ’dw'^W-V) 
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( W  ■  I ,  S,  Sj  . .  .  ) 
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3.2.3  Line  Loading 


The  plate  illustrated  in  Figure  12  is  considered 
next.  The  edge  conditions  are: 

x  =  0,a  - -  simply  supported, 

y  =  tb/2  .  »  clamped. 

The  deformation  equations  (89)  and  (90)  are  again 
used  to  give  the  forms  of  the  deflected  surface  for  the 
two  regions  of  y.  For  the  boundary  conditions  con¬ 
sidered  here,  the  constants  become: 
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~  XY  cosh 
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where  f*  is  defined  in  Paragraph  3*1*4  and  the 
parameters  T  and  <5  are  given  in  the  preceding 
paragraph . 

For  the  particular  case  of  the  line  load  F 
being  applied  along  the  x-axis, 

£  =0, 


and  the  constants  simplify  to: 


Mb 

L  oa  .  _  r 

z  (?-  V)  4 


3.3  Infinitely  Long  Plate 

An  infinitely  long  plate  with  the  long  sides 
simply  supported  is  now  considered.  Since  the  ratio 
of  b/a  is  very  large,  the  use  of  the  single  series, 
Levy-type,  solution  is  justified  and  the  deformation 
is  taken  in  the  form  of  Eq.  (49): 

**  I 
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and  Y.  satisfies  the  condition: 


Dry;-  2Dy 


o . 


(»»o) 


The  roots  of  the  characteristic  equation  have  four 
forms : 


roots  = 


+  fii  f^TW^ 

a  V  t>,  V  q1  Dr  ' 


Using  the  parameters, 


three  cases  (as  discussed  in  Paragraph  3.1.1)  must  be 
considered: 


Case  I  /A  ] 


Case  II  :  >U  -  i 


Case  III :  M  1 
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3.3.1  Line  Loading 


Let  the  line  load,  q,  act  along  the  x-axis  of 
an  infinitely  long  plate,  as  shown  in  Figure  18.  The 
deformation  will  be  taken  in  the  single  series  form, 
as  shown  by  Timoshenko  and  Woinowsky-Krieger  (Ref.  97). 


y 

Figure  18 

Infinite  Plate  -  Line  Load 


For  Case  I,  all  the  roots  of  Eq.  (ill)  are  real  and 
the  deformation, 
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For  Case  II,  there  are  two  double  roots,  and  the 
deformation, 


In  the  third  case,  the  deformation  has  the  form, 


The  following  notation  was  introduced  in  the  above 
expressions  for  the  sake  of  brevity: 


and  Q  is  the  load  per  unit  length. 


3.3.2  Concentrated  Force 


Now,  a  concentrated  load,  P,  is  applied  at 
(x  =  $  ,  y  =  0)  of  Figure  18.  Then,  the  deformation 
expressions  for  the  three  cases  of  elastic  properties 
become : 


Nowacki,  in  Reference  71,  obtains  a  solution  in 
closed  form  for  the  bending  mements  due  to  ft  load  of 
this  type. 

3.L  All  Edges  Clamped 

The  difficulty  presented  by  the  use  of  fixed 
supports  in  the  theory  of  plates  is  many  times  greater 
than  that  encountered  in  beam  theory.  The  usual  method 
of  solution  is  superposition.  Generally,  the  deformation 
is  found  for  a  simply  supported  plate  under  the  actual 
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loading,  and  it  is  added  to  the  deformation  due  to 
distributed  moments  along  the  edges. 

Any  of  the  previously  discussed  techniques  can  be 
applied  to  solving  this  problem  by  superposition. 

Since  Ando’s  method  has  been  used  rather  extensively  in 
this  chapter,  it  will  again  be  applied  here. 

The  problem  is  divided  into  two  parts:  1.)  Simple 
supports  at  x  =  ta/2  edges,  fixed  supports  at  y  =  tb/2 
edges,  and  under  actual  loading  (Fig.  19a);  2.)  Fixed 

supports  on  y  =  tb/2  edges  and  moment  distributed  along 
x  =  ta/2  edges  (Fig.  19b).  The  combined  effect  of  the 
two  parts  is  the  clamped  plate  shown  in  Figure  19c. 


Figure  19 

Clamped  Plate  by  Superposition 
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The  total  deformation,  w,  is  composed  of  two  terms 
corresponding  to  the  above  breakdown: 


w »  wa  ♦  ww 


where 


(>«) 


wa  =  deflection  of  plate  of  Figure  19a 
under  actual  loading, 

=  deflection  of  plate  of  Figure  19b. 

The  edge  moment  expression  is  taken  as: 

K  *  £  E„  >  <•"> 

where  ff*,  defines  its  intensity. 

From  this,  and  from  previous  considerations,  the 
deformations : 


wl  *  V  [Cc, «« ^  +  c4«*h  "®y) 

m 

4-  fc,  cosh  *t£E2lx  +  c*  cosh  '"I**)  sXy] , 

v  5  ^  <*  ^  *  b  j 


><«0 


91 


Evaluating  boundary  conditions,  the  total  deformation 
becomes : 


W  .  ♦  c,]  Cosh  -By 
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The  constants  of  Eq.  (122)  have  the  following  values: 
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The  constant,  A  ,  is  evaluated  for  each  particular 
loading  condition: 


a.)  Uniformly  distributed  load  (pQ): 


A 


4-P.c* 

D, n* 


(114) 


b.)  Concentrated  load  (P  -  at  distance  %  from  left 
y-edgej: 


2  Pa' 


jbXSl  . 


0*s) 


c.)  Hydrostatic  pressure  (maximum  p  -  varies  along 
x-axis ) : 


m-vt 

2(-0  ^mex 
£>m 


(>Z4) 


3 . 5  Summary 

Most  of  the  available  solutions  to  the  problem 
of  bending  of  rectangular  orthotropic  plates  require 
the  calculation  of  the  deformation,  w,  as  an  inter¬ 
mediate  step  prior  to  the  computation  of  desired 
stresses.  Regardless  of  the  method  used  to  arrive  at 
the  deformation  function,  the  expressions  given  in 
Part  1  are  applied  in  the  determination  of  stresses, 
forces,  and  moments. 

According  to  Eqs.  (27)  and  (23),  the  vertical 
shears , 
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The  moments,  from  Eqs.  (14),  (15),  and  (16)  are 
given  by: 
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M, 
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The  bending  and  shear  stresses,  respectively,  are 
given  by  Eqs.  (11),  (12),  and  (13)  as: 


2  a 

S<t 


The  meanings  of  the  parameters  used  in  the  above 
expressions,  and  the  coordinate  system  in  which  they 
are  oriented,  are  given  in  Part  1  of  this  report. 
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PART  4  -  ELASTIC  STABILITY 

A  great  deal  of  research  has  beer}  devoted  to  the 
problem  of  the  stability  of  various  types  of  stiffened 
plates,  with  reliance  upon  orthotropic  plate  theory. 

More  progress  has  been  made  in  buckling  theory  than  in 
bending  theory  for  two  reasons:  First,  the  buckling 
problem  of  aircraft  wings  has  been  of  greater  interest 
than  the  bending  problem,  and,  second,  the  solution  of 
the  former  is  considerably  simpler  than  that  of  trans¬ 
verse  bending. 

Some  of  the  authors  whose  work  is  of  importance 
due  to  either  original  contribution  or  summary  of 
available  methods  in  the  general  area  of  buckling  of 
orthotropic  plates  are  the  following:  Barg  (Ref.  5), 

Gerard  (Ref.  235),  Heck  and  Ebner  (Ref.  30),  Hu  (Ref.  36), 
Huber  (Ref.  37),  Kollbrunner  and  Meister  (Ref.  53), 
Lekhnitski  (Ref.  5&),  Nowacki  (Ref.  71),  Pfluger  (Ref.  73), 
Pochop  (Ref.  266),  Radok  (Ref.  267),  Sawczuk  (Ref.  $2), 
Schmit  and  Richer  (Ref.  317),  Schultz  (Ref.  S3),  Seide 
(Ref.  280),  Seydel  (Refs.  84  and  282),  Shuleshko  (Ref.  86), 
Sokolowski  (Ref.  284),  Strasser  (Ref.  285),  Thielemann 
(Ref.  95),  Timoshenko  (Ref.  289),  Timoshenko  and  Gere 
(Ref.  96),  Wilde  (Ref.  105),  Wittrick  (Ref.  107),  and 
Yusuff  (Refs.  Ill  and  300). 


preceding 
page  blank 
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Among  the  contributors  to  special  problems  are: 
Feldman  (Ref.  232),  who  presents  a  study  of  square 
orthotropic  plates;  Caldwell  (Ref.  214’),  who  considers 
the  special  case  of  a  plate  in  edge  compression  with 
reactions  provided  by  edge  shear;  Rockey  and  Cook 
(Refs.  269  and  270)  and  Symonds  (Ref.  287),  who  solve 
the  shear  buckling  problem;  Floor  and  Burgerhout 
(Ref.  233),  Hayashi  (Ref.  27),.  and  Wang  and  Zuckerberg 
(Ref.  294),  all  of  whom  investigate  the  post-buckling 
behavior  of  orthotropic  plates;  Klitchieff  (Ref.  250), 
who  limits  his  study  to  plates  with  longitudinal 
reinforcement . 

The  contributors  to  the  buckling  of  sandwich 
plates  are  not  mentioned  here.  A  summary  of  their 
work  is  contained  in  the  chapter  concerning  the  behavior 
of  sandwich  construction. 

This  chapter  contains  only  the  orthotropic 
buckling  theory.  It  should  be  noted  immediately  that 
this  theory  cannot  be  applied  indiscriminately  to  all 
stiffened  plates.  For  this  reason,  limitations  are 
discussed  first. 


4.1  Limitations 


Since  orthotropic  plate  theory  is  based  upon 

9 

average  flexural  and  twisting  rigidities  (Eqs.  (20), 
(21),  (22)),  it  can  be  stated,  in  general,  that  the 
theory  can  be  safely  applied  to  any  plate  with  a  large 
number  of  equal  and  equidistant  stiffeners  parallel  to 
one  of  the  edges  (Ref.  96). 

In  addition,  Gerard  and  Becker  (Ref.  236)  show 
by  test  results  that  orthotropic  theory  may  be  used 
in  the  following  cases:  a.)  plate  in  compression  with 
three  or  more  stiffeners;  b.)  plate  in  shear  with  any 
number  of  longitudinal  stiffeners;  c.)  any  plate  with 
transverse  stiffeners  which  has  a  low  or  high  value  of 

61 

bD 

where 

El  •  stiffness  of  stiffener, 
b  *  width  of  plate, 

D  •  rigidity  of  plate. 

4.2  General  Equation 

With  reference  to  the  derivation  of  Eq.  (31),  if 

the  transverse  load,  p,  is  now  replaced  by  in-plane 

forces,  N  ,  N  ,  and  N  ,  the  differential  equation  of 
x  y  xy 
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the  deformed  plate  is  written: 
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Or,  in  terms  of  the  stresses,  the  last  three  terms  of 
Eq.  (127)  can  be  written: 


+  °y 


(izft) 


where  t  is  the  thickness  of  the  orthotropic  plate, 
and  the  rigidities  are  defined  by  Eqs.  (20),  (21), 
(22),  and  (32). 

Buckling  due  to  the  compressive  loads  per  unit 
length,  N  and  N  ,  due  to  shear  per  unit  length,  N  , 

*  y 

and  due  to  combined  loadings  will  be  considered  here. 


Simple  Supports 

The  orthotropic  plate  with  simple  supports  along 
all  edges  (Fig.  20),  is  first  considered.  The  following 
paragraphs  give  solutions  for  various  loading  conditions. 


♦y  Figure  20 
Simply  Supported  Plate 


.1  Uniaxial  Compression 


Let  the  plate  of  Figure  20  be  loaded  uniaxially 
by  the  compressive  load  N  alone.  Then,  Eq.  (128) 

A 

becomes : 


The  deflected  surface  can  be  taken  in  the  form  of  a 


double  series, 
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But,  it  is  known  that  the  critical  condition  exists 
when  m  =  1,  n  =  1,  and  the  deformation, 

W  •  W.  S»M  H?  SIM  "7V  •  0*0 

a  b 

Substituting  this  into  Eq.  (129)  and  solving  for  the 
stress , 


<5 


b*t  L4  «l 


+  2D 
7 


0»2) 


The  smallest  value  for  the  critical  stress  is  obtained 
when 


and  the  critical  buckling  stress,  as  given  by  Timoshenko 
and  Gere  (Ref.  96), 

<C»  -  g'  [Joq  +  Dv]  OM) 

Lundquist  and  Stowell  (Ref.  253)  give  the  buckling 
stress,  for  the  condition  of  Eq.  (133),  as 

[2v^  D‘  +>4"  +  4<S|]  •  (BS) 
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Wittrick,  in  Reference  107,  gives  an  explicit 
solution  to  the  buckling  problem.  Using  Eq.  (130) 
as  the  equation  of  the  deformed  plate*,  he  takes  as 
the  solution, 


k  ■ 

L  m 

where 


(lit) 


k  »  buckling  coefficient, 

=  effective  aspect  ratio. 

The  buckling  coefficient  is  of  the  same  type  as  that 
used  for  isotropic  plates  and  can  be  applied  if  the 
effective  aspect  ratio, 

(?),  *  i 

is  employed. 

A  plot  of  Eq.  (136)  is  given  in  Figure  21.  From 
this,  the  coefficient,  k,  is  evaluated  for  a  given 
value  of  (a/b)e.  This  quantity  is  then  substituted 
into  the  expression  for  the  critical  buckling  load 
per  unit  length, 
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U, 3.2  Biaxial  Compression 


The  plate  of  Figure  20  is  now  loaded  by  compres¬ 
sive  loads  in  both  directions.  Thus,  Eq.  (127) 
becomes: 
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For  this  loading  condition,  Wittrick  (Ref.  107) 
applies  the  lower  curve  of  Figure  21,  but  the  effective 
aspect  ratio  must  be  modified  to  take  into  account  the 
load  Ny.  Thus,  the  effective  aspect  ratio  used  for 
this  case, 


i 


is  first  evaluated.  The  buckling  coefficient,  k,  is 
found  from  Figure  21,  and  the  following  expression  is 
used  to  complete  the  analysis: 


(Numbers  above  refer  to  order  of 
edge  supports  given  below. ) 
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Figure  21 

Buckling  Coefficients  for  Uniaxially  Loaded  Plates 

(from  Reference  107) 
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4.3.3  Edge  Shear 


When  the  plate  shown  in  Figure  20  is  acted  upon 
by  edge  shears  alone,  the  basic  equation  (127)  becomes: 


H  ^4w  +  2  D  +  ?  M  — w 

-  3k4  3x^yv  H  ZN»y  3x^y 
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Seydel  (Ref.  84)  also  makes  use  of  the  isotropic- 
type  buckling  coefficient  to  give  the  expression  for 
the  critical  shear  force  per  unit  length, 


(145) 


The  shear  buckling  coefficient,  k  ,  is  dependent  upon 

s 

the  effective  aspect  ratio  of  Eq.  (137)  and  upon  the 
parameter  D^/2  VD  xDy  .  The  coefficient  is  obtained 
from  Figure  22a  (which  is  given  by  Gerard  and  Becker 
in  Reference  236),  and  it  is  substituted  into  Eq.  (143). 

A  similar  result  is  given  by  Heck  and  Ebner  (Ref.  30) 


in  the  form: 


(144.) 


where  Cg  is  a  coefficient  dependent  on  the  parameters, 


« v  d, 


For  ranges 

0  «  E>i  1  .  0*1, 

the  curves  in  Figure  23  apply,  and  C  is  found  and 
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substituted  into  Eq.  (144). 
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Shear  Buckling  Coefficients 
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Figure  23 


Shear  Buckling  Parameter 


4.4  Two  Edges  Clamped.  Two  Edges  Simply  Supported 


The  rectangular  orthotropic  plate  with  two 
opposite  edges  clamped  and  the  other  two  edges 
simply  supported  (Figure  24)  is  now  considered. 


4.4.1  Clamped  Edges  Loaded;  Others  Unloaded 

When  N  acts,  and  N  =  N  =  0,  on  the  plate  of 
x  y  xy 

Figure  24,  Wittrick.  (Ref.  107)  gives  the  solution  as 
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where  k  is  the  coefficient  from  Figure  21  for  a 
particular  value  of  (a/b)e  given  by  Eq.  (137). 

4.4.2  Simply  Supported  Edges  Loaded;  Others  Unloaded 


When  N  acts,  and  N  =  N  =  0,  on  the  plate  of 
y  x  xy 

Figure  24,  the  solution  of  Wittrick  is  given  in  the  form 


047) 


and  k  is  available  from  Figure  21  for  the  effective 

aspect  ratio  given  by  Eq.  (137). 

Heck  and  Ebner  (Ref.  30)  give  as  the  solution  for 

the  case  of  D  >  D 
x  y 
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Lundquist  and  Stowell  (Ref.  253)  show  that  the 


critical  buckling  stress  for  this  case  is 


C^.l^U.SSAyqD,  +  l.217(M,D,+*(nD/)  ♦4.‘)43<3I] 
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for  the  minimum  condition, 
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4.4. 3  Biaxial  Compression 


When  N  and  N  act  on  the  plat  o  * gUre  24  and 
x  y 

N  =  0,  Wittrick  uses  the  curves  oj  r  ^ure  21,  but 
xy 

the  aspect  ratio  given  by  Eq.  (140)  must  be  applied. 
The  solution  is  obtained  from  Eq.  (141). 


4.4.4  Edge  Shear 

As  in  paragraph  4.3*3,  the  critical  shear  force 
per  unit  length  for  the  plate  of  Figure  24,  when 
N  =  N  =  0,  is  given  by 

*  y  _ 


where  the  shear  buckling  coefficient  is  available 
from  Figure  22b  (given  by  Timoshenko  and  Woinowsky- 
Krieger  in  Reference  (97)  for  plates  whose  edges  obey 
the  inequality, 


L.  5  All  Edges  Clamped 


A  plate  whose  four  edges  are  clamped  (Figure  25) 
is  considered. 
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Figure  25 
Clamped  Plate 


it, 5.1  Uniaxial  Compression 

When  N  is  acting  and  N  =  N  =  0,  the  critical 

x  y  xy 

load  on  the  plate  of  Figure  25  is  given  by  Wittrick 
(Ref.  107)  as 


(ISO) 
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and  the  buckling  coefficient  is  given  by  Figure  21 
for  the  effective  aspect  ratio, 

(2.)  =  2  A[Wi 
\b/€  b  \l  D.  ' 


4.5.2  Biaxial  Compression 

Timoshenko  and  Gere  {Ref.  96)  solve  the  problem 
of  a  biaxially  compressed  isotropic  plate  with  clamped 
edges,  assuming  the  deformation  function  as 

We  —  (|  -  Co  S  C°5  Z"T^)  *  (,5») 

and  using  strain  energy  to  calculate  the  critical 
buckling  stress.  A  similar  analysis  can  be  performed 
for  the  orthotropic  plate. 

4.6  Other  Buckling  Cases 

Several  other  cases  of  loadings  and  edge  con¬ 
ditions  are  considered  in  this  section. 


XV/ 
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.6.1  Combined  Shear  and  Compression 


Sandorff  (Ref.  315)  investigates  the  buckling  of 
a  plate  under  the  action  of  a  compressive  load  and 
shear  (Figure  26). 


ii 


Figure  26 
Combined  Loading 

The  assumption  is  made  that  b  is  infinite  and  an 
interaction  curve  (Figure  27)  is  presented.  The 
parameters  of  rhis  curve  are  stress  ratios.  That  i3, 
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Figure  27 

Interaction-*.  Shear  and  Compression 
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4.6.2  Loaded  Edges  Simply  Supported:  Other  Edges 
Clamped  and  Free 

An  implicit  solution  for  the  case  shown  in  Figure 
28  is  presented  by  Lundquist  and  Stowell  (Ref.  253) 
and  discussed  by  the  writer  (Ref.  309). 


The  critical  buckling  stress  is  given  by: 

2 

<£  =  (t j  D,  -  .Z70(^q  *  Mr,  D,)  4. 1 .712$] ,  OSS') 
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for  the  condition 


4.6.3  Three  Edges  Simply  Supported:  One  Edge  Free 

When  the  clamped  edge  of  Figure  28  is  replaced 
by  a  simply  supported  edge,  the  critical  buckling  stress 
of  Eq.  (154)  simplifies  to: 

for  the  same  limiting  condition. 
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PART  5  -  SANDWICH  PLATES 


The  name  "sandwich  plate"  describes  a  composite 
plate  which  consists  of  two  thin  sheets--faces--  and 
a  relatively  thick  core.  The  faces  are  generally 
of  aluminum  alloy  or  fiberglass,  while  the  core  can  be 
a  light  solid  material  such  as  plastic  and  wood  or  a 
thin  metal  sheet  formed  into  a  cellular  configuration 
("honeycomb").  The  basic  principle  of  sandwich  con¬ 
struction  is  the  reliance  upon  the  outer  faces  to 
carry  end  loads  and  upon  the  core  to  transmit  shear  and 
keep  the  faces  a  constant  distance  apart. 

Ordinarily,  the  modulus  of  elasticity  of  the  core 
material  is  extremely  small  (about  0.1%)  in  comparison 
to  that  of  the  faces.  Thus,  the  core  provides  very 
little  resistance  to  bending.  But,  because  its  modulus 
of  rigidity  is  also  very  low,  the  core  experiences 
appreciable  shearing  deformations.  Thus,  shearing 
deformations  must  be  considered  in  the  analysis.  In 
addition,  since  the  faces  or  the  core,  or  both,  may  have 
orthogonally  anisotropic  stretching  properties,  the 
composite  plate  may  well  be  orthogonally  anisotropic  in 
its  bending  properties.  Thus,  in  inany  cases,  the  theory 


1C  9 


which  should  be  used  to  analyze  sandwich  plates  is 
orthotropic  plate  theory,  with  consideration  cf 
shearing  deformations  (see  Chapter  1). 

The  amount  of  research  performed  in  the  field 
of  sandwich  construction  is  overwhelming.  For  this 
reason,  the  bibliography  of  this  report,  pertaining 
to  this  area,  is  by  no  means  complete.  For  specialized 
and  more  detailed  bibliographies  of  sandwich-plate 
literature,  see  References  131,  135,  and  137. 

Probably  the  most  important  of  the  papers  re¬ 
garding  the  theory  of  sandwich  plates  are  those  of 
Reissner  (Ref.  169), Hoff  (Refs.  139,  140),  and  Libove 
and  Batdorf  (Ref.  155).  Other  authors  who  have  made 
significant  contributions  in  the  investigation  of  the 
bending  and  the  buckling  of  sandwich  plates  are: 
Aleksandrova  (Ref.  113),  Anderson  (Ref.  116),  Anderson 
and  Updegraff  (Ref.  117),  Bijlaard  (Ref.  120),  Cheng 
(Ref.  125),  Dundrova  (Ref.  127),  Eringen  (Ref.  128), 
Goodier  (Ref.  132),  Guest  and  Solvey  (Ref.  136),  Hoff 
and  Mautner  (Ref.  142),  Hopkins  and  Pearson  (Ref.  143), 
Horvay  (Ref.  144),  Hubka,  Dow  and  Seide  (Ref.  145), 
Kimel  (Ref.  149),  Lewis  (Ref.  154),  Libove  and  Hubka 


(Ref.  156),  March  and  Smith  (Ref.  161),  Mushtari 
(Ref.  163),  Norris  (Ref.  164),  Raville  (Ref.  168), 

Seide  (Refs.  171,  and  172),  Seide  and  Stowell  (Ref. 

173),  Thurston  (Ref.  176),  Weikel  and  Kobayashi  (Ref. 
17&),  Williams  (Ref.  lSl),  Yen,  Salerno  and  Hoff 
(Ref.  1&4),  Yusuff  (Refs.  190,  191),  and  Zahn  and 
Cheng  (Ref.  192). 

Burns  (Ref.  123)  and  Burns  and  Skogh  (Refs.  212, 

213 )  present  minimum  weight  analyses  of  sandwich 
plates.  Kuenzi  (Ref.  152)  presents  design  criteria. 

Keer  and  Lazan  (Ref.  14S),  Mindlin  (Ref.  162),  and 
Yu  (Refs.  185-1S9)  consider  dynamic  behavior.  Manu- 
'  ture  of  sandwich  plates  is  discussed  by  Lowy  and 
Jaffee  (Ref.  157).  References  126  and  170  are  of 
considerable  value  to  the  designer  as  handbooks. 

This  chapter  is  intended  as  a  brief  summary  of 
some  of  the  highlights  of  the  small-  and  large-deflection 
theories  and  local  and  gross  instability  of  sandwich 
plates  from  several  of  the  papers  listed  in  Section  B 
of  the  Bibliography. 


5 . 1  Flexure 


The  analysis  of  sandwich  plates  in  bending 
follows  essentially  the  same  procedure  as  that  discussed 
in  Part  3*  However,  as  was  stated  previously,  the 
effect  of  transverse  shear  deformations  must  be  incor¬ 
porated  into  the  analysis. 

5.1.1  Small  Deflections 

Williams  (Ref.  181)  ,  Hopkins  and  Pearson  (Ref. 

143)  i  and  March  (Refs.  158-161)  were  among  the  first 
investigators  to  ac.junt  for  shear  deformations  by 
assuming  that  two  adjacent  vertical  strips  remain 
vertical  after  surface  displacement,  and  that  it  is 
the  shear  between  the  strips  which  limits  the  pene¬ 
tration  of  the  sinusoidal  face  displacement  of  the 
interface  of  core  and  facing.  Thus,  the  vertical  dis¬ 
placement  of  the  face  is  taken  in  the  form, 

V#  s  \  SIM  *2^  >  (I5S) 

while  the  displacement  of  the  core  in  the  same  direction 

is  given  by 
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The  constant,  k,  is  given  by 
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where 
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c  =  subscript  referring  to  core, 

E«  =  modulus  of  elasticity  in  z-direction, 

Gu=  shear  modulus  in  xz-plane, 

L  =  panel  length, 

w  =  number  of  half-waves  in  x-direction, 

*  =  longitudinal  coordinate, 

2  =  coordinate  through  thickness. 

More  general  and  rigorous  treatments  are  given 

by  Libove  and  Batdorf  (Ref.  155)  and  Reissner  (Ref. 

169).  Most  of  the  work  which  follows  chronologically 

is  based  on  the  principles  cited  in  these  papers. 

The  Libove-Batdorf  paper  (see  Paragraph  1.8) 

develops  a  theory  which  applies  to  orthotropic  or 

isotropic  plates  with  homogeneous  or  nonhomogeneous 

cores.  The  expressions  derived  therein  are  Eqs.  (45)  > 

(46)  and  (47).  Assuming  that  N  ,  N  ,  and  N  are 

x  y  xy 


constant  throughout  the  plate,  these  equations  are 
solved  for  the  deflection  and  the  shears: 


[  D]  w  e  -  [m]  f  » 
[d]  q.  •  -  [  N]  P  » 
[0]Qy  *-  t  p]  P  • 


(158) 


and  the  differential  operators  are  defined  in  the 
following  manner: 
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The  elastic  constants  which  appear  in  the  above 
equations  are  defined  by  Eqs .  (20),  (21),  (22),  (38), 
(39),  and  (40).  A  particular  problem,  thus,  is  solved 
by  using  Eqs.  (158)  for  a  given  set  of  boundary  con¬ 
ditions.  At  a  free  y-edge: 


Osn) 


M,  -  M,y  =  Q,  «  0. 


For  a  simply  supported  y-edge,  with  all  points 
along  boundary  prevented  from  moving  parallel  to 
edge : 

w  *  M  *  9  m  O  .  ( Ifco) 

For  a  simply  supported  y-edge,  with  all  points  free 
to  move  parallel  to  edge: 

w  *  =■  *  O  .  (ifoi) 

For  a  clamped  y-edge,  with  movement  parallel  to  edge 
prevented : 


w  -  (?w  _  ,  0 

KJ  \ 

For  a  clamped  y-edge,  with  all  points11'  free  to  move 
parallel  to  edge: 


Particular  results  for  the  bending  of  sandwich 
plates  under  uniform  and  concentrated  loads  are 
analyzed  by  Yen,  Gunturkun,  and  Pohle  (Ref.  183).  How¬ 
ever,  both  the  core  and  the  faces  are  considered 
isotropic  and,  therefore,  will  not  be  discussed  here. 


*Except  those  in  middle  surface. 
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The  same  is  true  of  Kuenzi*s  (Ref.  152)  report  on 
design  criteria. 

Raville  (Ref.  168)  presents  a  study  of  simply 
supported  rectangular  sandwich  plates,  with  isotropic 
facings  and  orthotropic  core,  acted  upon  by  a 
uniformly  distributed  load,  p.  The  expression  for 
the  maximum  deflection  (at  the  center  of  the  plate) 
is  given  by  a  double  Fourier  series: 


where 


s  = 


I  = 


vt,  v  z  ; 


(•64.) 


m,n  =  odd  integers, 

E  =  Young1 s  modulus  of  facings, 

M  =  Poisson’s  ratio  of  facings, 

9  =  a/b, 

while  a,b  and  the  coordinate  axes  are  consistent 
with  their  definitions  throughout  this  report. 

In  further  refinements  of'  small  deflection 
theory,  Hoff  (Ref.  139)  presents  a  variational 
approach  which  is  also  presented,  in  a  generalized 
fashion,  by  Eringen  (Ref.  128).  Both  papers  are 
limited  to  sandwiches  consisting  of  facings  and  core 
with  isotropic  properties. 

Cheng  (Ref.  125)  considers  the  problem  of  bend¬ 
ing  of  a  plate  with  isotropic  facings,  but  orthotropic 
core.  The  basic  equation  fb  r  bending  is  written: 


Ilf: 


in  which 
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h  = 
t  = 
E  = 

^  = 
A  = 


thickness  of  core, 
thickness  of  facings, 

Young1 s  modulus  of  facings, 
Poisson’s  ratio  of  facings, 
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The  homogeneous  solution  of  Eq.  (165)  satisfies 

[>-  -  d.  |;]aaw  =  o  ■ 


(«*) 
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This  form  is  used  to  solve  the  problem  of  a  simply 
supported  plate  under  a  uniformly  distributed  load 
(Raville’s  solution  is  given  elsewhere  in  this  section). 
Since  the  load  intensity,  p,  is  constant,  the  right 
hand  side  of  Eq.  (165)  reduces  to  p/D,  and  its  parti¬ 
cular  solution  is  taken  in  the  form: 


•s 


The  maximum  deflection,  at  the  center  of  the  plate, 
is  then  given  by: 
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5*1.2  Large  Deflections 


Generally  speaking,  the  study  of  large  deflections 
of  sandwich  plates  is  an  extension,  through  the 
inclusion  of  shear  effects,  of  von  Karmanfs  nonlinear 
theory  of  elasticity.  Such  a  study  is  quite  important 
because,  since  the  advent  of  high  strength  alloy  steels 
and  titanium  alloys,  it  is  possible  to  construct  a 
sandwich  consisting  of  a  soft  core  and  extremely  thin 
facings . 

Reissner  (Ref.  169)  considers  a  plate  with 
isotropic  facings  and  core,  and  uses  the  simplifying 
assumption  that  stresses  in  the  core,  whose  vectors  lie 
in  the  xy-plane  (see  Fig.  1),  and  the  variation  of 
face  stresses  over  the  thickness  of  the  faces,  can  both 
be  neglected.  In  other  words,  the  face  sheets  act 
as  membranes,  while  the  core  resists  transverse  shear 
and  normal  stress.  Introducing  an  Airy  stress  function, 
F,  two  simultaneous  differential  equations  are  written 
to  describe  the  behavior  of  an  isotropic  sandwich  plate: 


aaF 


04-0 
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The  nonlinear  theory  of  a  symmetrically  loaded 
sandwich  beam  is  considered  by  Zahm  and  Cheng  in 
Reference  192,  based  on  the  assumptions  that  the 
exact  load  distribution  is  known  and  can  be  repre¬ 
sented  by  a  Fourier  series  and  that  the  core  material 
is  anti-plane. 

5.2  Elastic  Stability 

Buckling  of  sandwich  plates  is  divided  here  into 
two  parts,  local  instability  and  general  (gross) 
buckling. 

5.2.1  Local  Buckling 

5. 2. 1.1  Face  Wrinkling 

Several  investigators  consider  local  crinkling 
of  the  faces  in  a  plate  consisting  of  isotropic 
facings  and  an  orthotropic  core.  Among  those  who 
originated  such  studj.es  are  Williams  (Ref.  181)  , 

Hoff  .and  Mautner  (Refs.  141,  142),  and  Goodier  (Ref. 
132).  A  relatively  simple  formula  for  the  critical 
compressive  load  (applied  in  x-direction)  per  unit 


width  of  panel  is  given  by  Williams: 


where  the  subscript,  c,  denotes  ’’core",  the  coordinates 
are  those  of  Figure  1,  and  all  other  symbols  are  as 
defined  in  Part  5»1. 

Goodier  and  Hsu  (Ref.  133)  investigate  face 
wrinkling  in  a  non-sinusoidal  mode  and  derive  a  lower 
critical  load.  Yusuff  (Refs.  190,  191)  considers 
initial  imperfections  and,  for  the  plate  shown  in 
Figure  29,  derives  equations  for  critical  stress,  based 
on  local  wrinkling,  of  a  plate  with  isotropic  faces 
and  core. 


Figure  29 

Wrinkling  of  Facing 
W  is  the  "zone  of  displacement", 


For  a  failure  of  the  core  by  tension  or  compression, 
the  critical  stress, 


(I7») 


For  core  failure  due  to  shear 


The  ratio,  AQ/F,  is  known  as  the  "waviness  parameter 
and  is  given  in  Figure  30.  The  curve  given  in  this 
figure  has  been  confirmed  by  various  tests,  among 
them  those  of  Harris  and  Crisman  (Ref.  138). 
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Figure  30 

Waviness  Parameter  versus  Core  Moduli 

5 >2.1.2  Intracellular  Buckling 

A  second  form  of  local  instability  is  that  of 
"dimpling”  or  intracellular  buckling.  This  failure,  of 
course,  is  dependent  upon'  the  core  configuration. 

Anderson  and  Updegraff  (Ref.  117)  investigated 
the  buckling  of  truss-core  sandwich  panels  experimentally. 
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where  GQ  is  the  shear  modulus  of  the  core  and  the 
definitions  of  all  other  terms  are  given  in  Eq. 
(165). 

These  equations  (127  and  123),  together  with  the 
boundary  conditions,  determine  the  functions  F  and  w. 
For  the  solution  of  such  differential  equations,  the 
reader  should  refer  to  page  425  of  Reference  96. 

The  large  deflections  of  a  sandwich  plate  with 
orthotropic  core  between  isotropic  face  sheets  are 
considered  by  Alwan  (Refs.  114,  115).  Here,  the 
solution  is  given  by  two  differential  equations,  the 
first  of  which  is  Eq.  (169)  and  the  second, 
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Writing  the  critical  buckling  stress  of  the  faces  in 

i  ►  '  .  ■ 

the  intracellular  mode, 
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where  the  symbols  are  defined  in  Figure  31  and  in 
Part  5,1.  The  values  of  the  buckling  coefficient,  k, 
are  given  in  Figure  32. 
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Figure  31 

Truss-Core  Sandwich 
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Figure  32 

Truss-Core  Buckling  Coefficient 


A  design  curve,  for  honecombs,  based  on  intra¬ 
cellular  buckling,  is  given  on  page  C12.8  of  Reference 
122.  When  uniaxial  compression  is  applied  in  the 
direction  of  the  ribbon, 


where 


*1 

n 


is  the  plasticity  reduction  factor, 


.  -iix- 
E  +  Et 


and  E,p  is  the  tangent  modulus  of  the  facings.  The 
geometry  of  the  honeycomb  is  given  in  Figure  33  • 
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Figure  33 

Honeycomb  ] ] 

T| 

5. 2. 1.3  Shear  Crimping  -* 

The  third  form  of  local  instability  is  excessive  .1 

transverse  shear  stress  (shear  crimping).  Such  a  -j 

failure  may  occur  when  the  sandwich  panel  is  not 
initially  perfectly  flat.  )) 
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The  problem  Is  considered  by  Boiler  and  Norris 
(Ref.  121)  and  Bruhn  (Ref.  122).  The  latter  gives 
the  design  equation, 


07S> 


i 

with  terms  defined  in  the  same  manner  as  those  of 

Eq.  (177) .  j 

5.2.2  General  Buckling 

The  general  buckling  of  sandwich  plates  has  been 
investigated  by  several  authors  for  many  types  of  end 
conditions  and  loadings. 

\ 

t 

5. 2. 2.1  Simply  Supported  Edges 

Buckling  of  simply  supported  sandwich  plates  is 
discussed  by  Gerard  and  Becker  (Ref.  236),  Hoff 
(Refs.  139»  140),  Kimel  (Ref.  149 ),  March  and  Smith 
(Ref.  161),  and  Seide  and  Stowell  (Ref.  173). 
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For  the  isotropic  sandwich  plate  in  uniaxial 
compression ,  the  formula  for  the  critical  buckling 
stress  is  given  by  Hoff: 


QjJ  is  the  buckling  stress  of  the  unconnected 


faces, 


Hit' 


> 


F  is  the  form  factor, 

F  -  l  +  3(l+£)1  . 

and  C  is  a  reduction  factor  for  the  buckling  of  the 
composite  panel.  Its  values  are  given  by  the  curve 
in  Figure  34»  for  various  parameters, 


(m) 


R 


Reduction  Factor  -  Simple  Supports 

Kimel,  in  Reference  149,  performs  a  theoretical 
analysis  of  the  simply  supported  plate  under  the 
combination  of  edgewise  bending  and  compression  on 
opposite  edges.  The  sandwich  considered  is  made  up 
of  isotropic  facings  and  an  orthotropic  core.  The 
results  of  this  analysis  are  quite  lengthy  and  will  not 
be  presented  here. 
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A  uniaxially  compressed  sandwich  with  an 
orthotropic  core  is  discussed  by  March  and  Smith 
(Ref.  l6l).  With  the  load  applied  along  the  a-edges, 
the  critical  buckling  load, 


06©) 


where 
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q 
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E, 


n  =  number  of  half-waves  in  x-direction. 


1  TM 


A  family  of  curves,  for  various  values  of  n,  must 
be  constructed  in  order  to  make  Eq.  (180)  suitable 
for  design. 

5. 2. 2. 2  Loaded  Edges  Pinned.  Other  Clamped 

This  case  is  considered  by  Yen,  Salerno,  and 
Hoff  (Ref.  184).  A  close  approximation  to  the  true 
value  of  the  buckling  load  is  obtained  by  establishing 
a  lower  bound  by  Leggett’s  method  and  ^n  upper  bound 
by  Galerkin’s  approach.  The  derived  expressions  are 
extremely  complex,  but  a  diagram,  well  suited  for 
design,  is  presented. 

A  Metalite  sandwich  (isotropic  facings,  ortho¬ 
tropic  core)  is  studied  by  Seide  (Ref.  171).  In  this 
paper,  the  flexural  rigidity  of  the  faces  is  neglected, 
thus  underestimating  the  shear  stiffness  of  the  core 
(discussed  by  Bijlaard  in  Reference  120).  Using  a 
combined  result  of  Seide  and  Bijlaard  (modification 
of  former’s  method  by  including  face  rigidities),  the 
critical  buckling  load  along  the  pinned  edges, 
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where 

K. 

= 

buckling  coefficient  from  Figure  35, 
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= 

plasticity  factor  (see  5. 2. 1.2), 
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and  all  other 

terms  are  those  used  with  consistency 

throughout  this  paper. 

March  and  Smith  (Ref.  161)  take  the  x-edges  as 
loaded  and  simply  supported  and  derive  the  expression 
for  the  critical  buckling  load  of  a  sandwich  with  an 
orthotropic  core  as: 


(\%l) 


where  all  terms  are  defined  in  Eq.  (180). 

5. 2. 2. 3  Loaded  Edges  Clamped.  Others  Pinned 

\ 

The  solution  to  the  problem  of  the  x-edges 
loaded  and  clamped  and  the  y-edges  simply  supported  is 
given  by  March  and  Smith  (Ref.  161)  for  the  sandwich 
with  an  orthotropic  core : 
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where,  for  one  half-wave, 
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for  two  half-waves, 
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for  three  half-waves, 


for  four  half-waves, 
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♦  Mill  , 


and  all  terms  are  defined  in  Eq.  (180). 


5. 2. 2. 4  All  Edges  Clamped 

March  and  Smith  (Ref.  161)  solve  the  problem 
of  the  clamped  sandwich  plate  with  isotropic  facings 
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and  orthotropic  core.  The  critical  buckling  load 
along  x-edge, 
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in  which  the  buckling  coefficient  for  one  half-wave, 


for  two  half-waves, 
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for  three  half-waves, 
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for  four  half-waves, 
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and  all  terms  are  defined  in  Eq.  (180). 
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Buckling  Coefficients  for  Eq.  (18 
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This  problem  is  also  solved  by  Thurston 
(Ref.  176).  The  buckling  load  is  given  by: 


~irxEt  (v.+f)1, 
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in  which  K*  is  a  buckling  coefficient  given  by 
Figure  36.  In  the  Figure,  S  is  the  sandwich 
parameter, 


S 


(iM) 


where 

D  .  et  (hitf  . 
Z  (1-aO 


111 


Figure  36 

Buckling  Coefficients  for  Eq.  (185) 

5.3  Summary 

Only  the  most  common  cases  of  loading  and  end 
supports  have  been  considered  in  this  chapter.  For 
other  combinations  of  end  conditions  and  loads,  the 
reader  is  advised  to  search  the  titles  of  Section  B 
of  the  Bibliography. 
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PART  6  -  CORRUGATED  PLATES 


Corrugated  plates  were  first  introduced  as  a  node 
of  construction  in  England  in  the  late  1920* s.  Their 
first  use  was  in  the  form  of  tank  and  sub-division 
bulkheads  in  ships.  Corrugated  sheet  steel  was  later 
used  in  the  construction  of  aircraft  wings  and  fuse¬ 
lages.  In  recent  years,  corrugated  sheet  attached  to 
a  flat  plate  (generally  called  "corrugation-stiffened 
panel")  has  been  used  in  the  forming  of  skin  of  high¬ 
speed  flight  vehicles  such  as  the  X-15.  A  more  recent 
advance  has  been  the  adoption  of  corrugated  plate  for 
use  as  core  in  an  all-metal  sandwich. 

6.1  Flexure 

The  flexure  of  corrugated  plates  of  various 
configurations  has  been  discussed  by  several  authors. 
Among  these,  the  following  are  the  most  prominent: 
Andreeva  (Ref.  193),  Caldwell  (Refs.  196,  197), 

Shibuya  (Ref.  205),  Stroud  (Ref.  335),  Platus  and 
Uchiyama  (Ref.  204),  Bergman  and  Reissner  (Ref.  195), 
and  Timoshenko  and  Woinowsky-Krieger  (Ref.  97). 

The  discussion  which  follows  is  divided  into  two 
parts.  First,  a  single  corrugated  plate  is  considered, 
and,  second,  a  built-up  configuration  of  corrugated 
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and  flat  plates  Is  discussed. 


6.1.1  Single  Sheets 


The  bending  of  single  corrugated  sheet  is  taken 
up  by  Timoshenko  and  Woinowsky-Krieger  (Ref.  97). 

For  a  sinusoidal  configuration  (Fig.  37),  given  by 
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the  orthotropic  constants  are: 
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E  =  modulus  of  elasticity  of  material, 
M  =  Poisson’s  ratio  of  material, 
s  =  arc  length  of  half-wave. 


Figure  37 

Sinusoidal  Corrugation 


From  the  geometry,  the  arc  length  and  the  moment  of 
inertia  are  approximately  given  by 


The  orthotropic  plate  theory,  discussed  in  Part  3,  is 
than  applied  for  the  solution  of  a  particular  bending 
problem. 

Corrugations  of  a  different  geometry  (Fig.  38) 
are  considered  by  Caldwell  (Ref.  196).  The  author 
takes  into  account  the  effect  of  inelastic  behavior. 


Figure  38 

’’Hat -type"  Corrugation 
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Experiments  show  that  elastic  stresses  and  deformations 
can  be  calculated  from  ordinary  bending  theory  using 
a  moment  of  inertia, 

1  -  +  O  > 

<o 

and  a  section  modulus, 

2  -  (it>4c)  .  0*4) 

3 

Assuming  that  the  strains  are  large,  the  ultimate 
bending  strength  of  the  panel  is  given  by: 


where  is  the  yield  strength  of  the  material. 


6 . 1 .2  Built -up  Panels 


Stroud  (Ref.  335)  calculates  the  elastic  con¬ 
stants  of  a  corrugated  plate  attached  on  one  side 
to  a  flat  plate  (Fig.  39). 
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Figure  39 

Corrugation -Stiffened  Plate 


The  differential  equation  of  this  orthotropic  plate 
is  written  in  a  form  slightly  different  from  that  of 
Eq.  (31): 
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is  the  bending  stiffness  in  a  direction  parallel 
to  the  troughs  (’’into  the  paper”  -  -  -  Fig..  39), 


and,  with  E  being  the  modulus  of  elasticity  of  the 
material,  its  value  is: 
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in  which  the  location  of  the  neutral  axis  is  given  by 
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The  bending  stiffness  in  a  direction  normal  to  the 
trough  axis , 


Dr  - 


a  4  U 

q-  q 


0181 


148 


«  » 


where 


_ l_ 

i  ♦  (  R,^  ©  4  k  ♦  ^ 


*“ H  [R*  (9- *••*€?)  4  jk4-  £**im  ©  4*  .  • . 

» •  • 

4  Rt  (nfi  -  R^©  coj  G  4  <?,>•»!  ©)  4*  J  43 

**£*>]  ('-'V*5  ' 

-  <!-*»  **  )  , 


l?V  >»*9)  t  j>c  4l^Q4  ... 

R,*  (*£  -  2  ft*  ©  4  ft*  if  )  4  4  .  .  . 

»  Ki  (w9  -  (?tOc<40  4  I?,,  im  e1)  4.  ... 

♦  j kv4»a  G>  ■*  y  *«*'•  4  (2,1^0  -  ... 


2n  Ru  3  Co*  ©  4  2  **  ^  iM  G  4  •  •  • 


The  twisting  stiffness,  with  G  being  the  shear  modulus 
of  the  material, 
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The  bending  moments  are  evaluated  by  substituting  the 
above  values  (in  which  /4  is  the  Poisson’s  ratio  of  the 
material)  into 
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The  same  problem  is  considered  by  Wempner  and 
McKinley  (Ref.  130)  and  similar  results  are  derived 
for  plates  connected  by  one  weld-line  (a  =  0  in  Fig.  39). 


6.2  Elastic  Stability 


The  buckling  of  corrugated  plates  has  been  studied 
by  various  investigators.  Among  these  are  Andreeva 
(Ref.  193),  Ashwell  (Ref.  194),  Libove  and  Hubka 
(Ref.  156),  Plaineveaux  (Ref.  203),  Stroud  (Ref.  326), 
Seydel  (Ref.  84),  Timoshenko  and  Gere  (Ref.  96),  Wempner 


and  McKinley  (Ref.  180),  and  Dean  (Ref.  198). 

This  portion  will  again  be  divided  into  two 
phases:  the  first  dealing  with  single  corrugated 
sheets  and  the  second  with  panels  consisting  of  flat 
plates  attached  to  corrugated  sheets. 

6.2.1  Single  Sheets 

Since  consideration  of  bending  requirements  alone 
would  dictate  the  design  of  a  *hin  plate  with  deep 
troughs,  it  is  necessary  to  consider  limitations  im¬ 
posed  upon  such  a  design  by  local  instability  of  the 
component  "walls”.  Caldwell  (Ref.  196)  discusses  this 
problem,  taking  into  account  the  interaction  between 
webs  and  flanges  and  referring  all  symbols  to  the  con¬ 
figuration  of  Figure  38. 

/ 

The  critical  buckling  stress ,  based  on  local 
instability  (and,  thus,  independent  of  panel  length 
and  edge  supports),  is  given  by: 


Q  p 

where  D  is  a  material  property,  Et;/12(l -)U  ) ,  and  Kc 
is  a  constant  dependent  upon  the  dimensions  of  the 
corrugation  (see  Figure  38  and  40). 
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Timoshenko  and  Gere  (Ref.  96)  discuss  the  work 
of  Bergmann  and  Reissner  (Ref.  195)  and  Seydel  (Ref.  84) 
for  the  critical  buckling  shear  load,  and  give  the 
expression : 

(Mj  .  4C,  £5$  ,  <»•>> 

where  C  is  the  buckling  coefficient  given  in  Figure  23. 
The  equation  (202)  is  that  given  in  paragraph  4*3.3  and 
all  definitions  are  given  there. 


Figure  40 

Buckling  Coefficient  of  Eq.  (201) 
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6.2,2  Built-up  Panels 


The  formulas  derived  by  Stroud  (Ref.  335)  can  be 
extended  to  the  buckling  case  be  rewriting  the  differ¬ 
ential  equation  (196)  such  that  it  includes  edge 
loadings : 
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Thus,  if  the  orthotropic  constants  defined  by  Eqs. 
(197),  (19S),  and  (199)  are  substituted  into  Eq.  (203), 
the  solutions  given  in  Part  4  can  be  applied  to  par¬ 
ticular  problems. 

This  problem  is  also  discussed  by  Wempner  and 
McKinley  (Ref.  ISO)  and  theoretical  results  are  com¬ 
pared  with  experimental  data. 


6.3  Summary 


The  increase  in  load-carrying  ability  accomplished 
by  constructing  corrugations  parallel  to  the  loading 
and  the  high  flexural  rigidity  of  corrugated  plates  are 
well  recognized.  However,  plates  of  this  type  have 
several  disadvantages:  a)  i.e.  strength  is  increased 


in  one  direction,  while  it  is  decreased  in  the  direction 
normal  to  it;  b)  weight  per  unit  projected  area, is 
greater  than  that  for  flat  plates;  c)  workability  of 
corrugated  sheet  is  limited. 

For  these  reasons,  various  investigators  have 
explored  other  configurations  which  can  be  formed  as 
easily  as  corrugations.  Smith  and  Gray  (Ref.  320) 
discuss  a  method  for  redistributing  metal  in  flat 
sheets  for  the  purpose  of  obtaining  maximum  strength- 
weight  ratios  and  rigidity -weight  ratios.  Farmer  and 
Spangler  (Ref.  308)  investigate  a  pattern  of  "calottes”, 
or  hemispheres,  pressed  into  a  flat  sheet  in  a  random 
pattern.  This  writer  (Ref.  309)  and  Tewes  (Ref.  321) 
present  a  study  of  plates  stiffened  by  "dimples" 
(L-shaped  depressions  pressed  into  the  sheet  in  an 
interlocking  manner).  "Waffle  beading",  forming  of 
intersecting  sinusoidally  shaped  beads  in  a  sheet 
metal  panel,  is  discussed  by  Rogge  in  Reference  314. 


PART  7  -  RIB -REINFORCED  PLATES 


The  bending  and  buckling  of  flat  plates  stiffened 
by  longitudinal  or  transverse  ribs,  or  both,  is  a 
subject  which  has  received  much  attention  in  the 
literature.  Early  studies  were  concerred  with  multi¬ 
piece  skin-and-stiffener  combinations,  while  more 
recent  investigations  are  generally  confined  to  in¬ 
tegrally-stiffened  structures.  Both  types  are  con¬ 
sidered  in  this  chapter. 

7.1  Flexure 

Among  the  many  investigators  concerned  with  the 
bending  of  rib-stiffened  plates  are  the  following: 
Borcz  (Ref.  326),  Chapman  and  Slatford  (Ref.  217), 
Floor  and  Burgerhout  (Ref.  233),  Giencke  (Refs.  237- 
239),  Giuliani  (Ref.  240),  Holmes  (Ref.  244),  Huber 
(Ref.  40),  Jaeger  and  Hendry  (Ref.  246),  Lekhnitsky 
(Ref.  5#) ,  Richart,  Newmark  and  Siess  (Ref.  268), 
Schade  (Refs.  273-276),  Schumann  (Refs.  278,  279), 
Terazawa  and  Tarada  (Ref.  288),  Timoshenko  and 
Woinowsky-Krieger  (Ref.  97),  Tolotti  and  Grioli 
(Ref.  290),  and  Yaraaki  (Ref,  29S). 


I 


The  bending  of  orthogonally  stiffened  plates 
(to  which  this  study  is  limited)  can  be  analyzed  by 
means  of  orthotropic  plate  theory,  as  described  in 
Parts  1  and  2.  The  elastic  constants  of  the  govern¬ 
ing  differential  equation  (31)  must  be  determined 
experimentally,  according  to  the  methods  outlined  in 
the  following  chapter.  Then,  deformations,  stresses, 
etc.,  can  be  evaluated.  Orthotropic  constants  for 
several  cases  of  stiffening  are  discussed  here. 

Lekhnitsky  (Ref.  53)  considers  a  plate  with 
equally  spaced  one-directional  stiffeners  symmetrical 
with  respect  to  the  midplane  of  the  plate  (Figure  41). 
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The  constants  of  the  characteristic  equation  (31)  are: 

D  -  D  •  -  , 

H  * i  12  0-x1-) 

D  -  +  E±±Ll)'  . 

1  iz(i-,u«)  L,  -* 

where  subscripts  p  and  s  refer  to  properties  of  plate 

and  stiffener,  respectively. 

If  the  plate  of  Figure  41  is  now  reinforced  by 

another  set  of  equal  stiffeners  in  the  y-direction, 

a  distance  apart, 

d  .  liil  *  E-iiL)« , 

»  L, 

D  «  +  J  * 

1  L. 

D  -  -Ir*-  , 


<2eS) 
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where 


Ef>E,  m  moduli  of  elasticity  of  plate  and 
stiffeners,  respectively, 


M  •  Poissons  ratio  of  plate  material, 

(i )  *  moment  of  inertia  of  one  stiffener, 
**  parallel  to  x-axis, 


(I  )  ■  moment  of  inertia  of  one  stiffener, 
*  parallel  to  y-axis. 
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The  constants  for  a  plate  integrally  reinforced 
by  equal  and  equidistant  ribs  are  derived  by  Pfliigger 
(Ref.  264)  and  Timoshenko  and  Woinowsky-Krieger 
(Ref.  97).  The  latter  give  the  following  expressions 
for  the  plate  shown  in  Figure  42: 


Integrally 


Figure  42 

Stiffened  Plates  Governed  by  Eqs. 


(206) 


where 

I,  *  moment  of  inertia  of  one  T-sectidn 
of  width  L  , 

C  ■  torsional  rigidity  of  one  rib, 

CjJ  *  torsional  rigidity  of  plate  without 
ribs , 

OC  -  . 

In  addition  to  solutions  which  make  use  of 
orthotropic  theory,  in  which  the  interaction  between 
plate  and  stiffener  is  generally  neglected,  several 
writers  presented  more  sophisticated  direct  methods  of 
attack. 

One  of  the  early  investigators,  Schade  (Ref.  276), 
considers  a  rectangular  plate  with  longitudinal  and 
transverse  stiffeners,  and  with  the  two  y-edges  clamped 
and  the  x-edges  free,  acted  upon  by  a  uniformly  distri¬ 
buted  load,  p  (Figure  43). 
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The  deformation  expression,  taken  as  a  single 
series,  has  the  form: 


of  longitudinal  stiffener, 


•  • 

A'l  *  Xfn  *  ***  «  average  unit  moment  of  inertia 

2» 

of  plate  alone, 

I*  85  moment  of  inertia  of  longitudinal 
stiffener  and  effective  skin, 

ly  »  moment  of  inertia  of  transverse 
stiffener  and  effective  skin, 

*  moment  of  inertia  of  central 
longitudinal  stiffener  and  effective 
skin, 

Iy 

*1  -  -g  -  unit  moment  of  inertia  of  plate 
alone. 


Figure  43 

Stiffened  Plate  Governed  by  Eq.  (207) 


A  numerical  procedure  for  the  direct  calculation 
of  stresses  is  presented  by  Kempner  (Ref.  248).  Holmes 
(Ref.  244)  discusses  the  effect  of  shear  lag  upon  the 
bending  of  stiffened  plates. 
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7.2  Elastic  Stability 


> 


The  problem  of  buckling  of  stiffened  plates  is 
extremely  important  to  the  designers  of  aircraft 
because  of  the  rigid  requirement  of  aerodynamic 
smoothness  of  wing  surfaces.  For  this  reason,  this 
area  has  been  investigated  by  a  great  number  of 
researchers.  Among  the  more  important  contributors 
are  the  following:  Burns  and  Skogh  (Ref.  212), 

Caldwell  (Ref.  214),  Dow,  Hickman  and  Rosen  (Ref.  228), 
Gerard  (Refs.  234,  235),  Gerard  and  Becker  (Ref.  236), 
Irving  and  Mullineux  (Ref.  245),  Klitchieff  (Ref.  250), 
Nowacki  (Ref.  262),  Pfluger  (Ref.  264),  Pochop  (Ref.  266), 
Radok  (Ref.  ,267),  Rockey  and  Cook  (Refs.  269,  270), 
Rosenhaupt  (Ref.  271),  Seide  (Ref.  280),  Seydel  (Ref.  282), 
Sokolowski  (Ref.  284),  Strasser  (Ref.  285),  Timoshenko 
(Ref.  289),  Timoshenko  and  Gere  (Ref.  96),  Tsuiji 
(Ref.  291),  Wang  and  Zuckerberg  (Ref.  294),  and  Yusuff 
(Ref.  300). 


When  orthotropic  theory  applies ,  the  deflections 
of  the  stiffened  plate  are  governed  by  Eq.  (127): 
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and  the  values  of  the  constants  Dx,  D^,  and  are 
given  by  Eqs.  (204)  or  (205)  or  (206)  or  can  be  de¬ 
termined  by  experimental  means,  using  methods  outlined 
in  the  next  part.  The  expressions  derived  in  Part  4 
can  then  be  used  to  solve  particular  problems. 

A  very  high  percentage  of  papers  dealing  with  this 
area  are  experimental  in  nature,  and,  thus,  solutions 
are  generally  given  in  the  forms  of  curves.  The  ex¬ 
tremely  large  amount  of  information  available  renders 
impossible  a  summary  of  basic  expressions  and  data  in 
this  report.  Since  the  literature  contains  data  on 
many  different  configurations  of  stiffened  panels,  it 
is  impossible  to  present  a  concise  guide  here. 

For  this  reason,  the  analyst  is  referred  to  the 
bibliography  of  this  paper.  The  titles  of  most  papers 
adequately  describe  the  specialized  topics  with  which 
they  deal.  For  excellent  summaries,  Handbook  of 
Structural  Stability.  Parts  V  and  VII  (Refs.  235  and 
236),  are  particularly  recommended.  Bruhn’s  well- 
known  reference  (Ref.  122)  also  contains  a  great  deal 
of  information  dealing  with  the  buckling  of  stiffened 
panels . 
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7.3  Summary 

The  flat  plate  reinforced  by  either  longitudinal 
or  transverse  stringers,  or  both,  remains  an  extremely 
important  problem  in  the  design  of  ships,  flight 
vehicles,  and  bridges.  Only  the  static  behavior  has 
been  discussed  here.  However,  there  are  several  in¬ 
vestigators  who  have  also  concerned  themselves  with  the 
dynamic  problem.  Among  these  are:  Feldman  (Ref.  232), 
Greenspon  (Ref.  23),  and  McElman,  Mikulas  and  Stein 
(Ref.  256). 
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PART  8  -  DETERMINATION  OF  ELASTIC  CONSTANTS 


In  the  previous  chapters,  use  is  made  of  ortho¬ 
tropic  plate  theory  to  solve  various  bending  and 
buckling  problems.  In  order  that  this  theory  may  be 
applied,  it  is  first  necessary  to  determine  the  various 
elastic  constants  of  orthotropy.  Let  us  review  the 
characteristic  differential  equation  for  a  general 
loading: 
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In  this  expression,  the  shear  deformation  effects  are 
neglected.  Since  such  a  simplification  is  permissible 
in  many  cases,  the  determination  of  the  three  constants, 
Dx,  Dy,  and  D^,  will  first  be  considered. 

In  the  more  refined  case,  where  shear  deformations 
are  accounted  for  (see  Paragraph  1.8)  three  additional 
constants  must  be  evaluated.  The  methods  of  evaluation 
of  all  of  these  values  are  discussed  in  this  part. 

8.1  Bending  Stiffness 

The  original  scheme  for  measuring  the  bending 
stiffness  of  orthotropic  plates  is  given  by 


4 


Bergstrasser  (Ref.  325).  The  same  method  is  then 
discussed  by  Thielemann  in  1945 ,  and  later  by  Hearmon 
and  Adams  (Ref.  29),  Hoppmann  et  al  (Refs.  329-332), 
and  Tsai  and  Springer  (Ref.  338). 

First,  the  relationship  between  stresses  and 
strains,  Eq.  (3),  is  recalled: 
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Since  S12  =  S21,  there  are  four  elastic  constants  to 
be  determined  in  the  case  where  shearing  deformations 
are  neglected. 

For  the  determination  of  the  bending  constants : 
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plate  is  oriented 
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and  supported  as  shown  in  Figure 


44  • 
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Figure  44 

Test  for  Bending  Stiffness  (I) 

Thus,  a  uniformly  distributed  bending  moment,  Mx,  is 
applied  to  the  plate.  Therefore,  the  deformation  at 
any  point  (x,  y), 

w  -  +  sltyl)  +  Ax  +  By  +  C  >  Czos) 

where  A,  B,  and  C  are  constants.  The  moment  in  terms 
of  the  total  applied  load,  P, 
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(20*0 


From  this,  the  deflection, 

w*  oP  (l  “  f  )(SHXl  +  Sityt)  +  +  By  +  C  .  (210^ 


It  is  now  necessary  to  make  two  deflection  measurements 
in  order  to  evaluate  the  constants  and  S^.  These 
two  measurements  are  at 

x  =  0,  y  =  0, 
x  =  0 ,  y  =  ^ . 

It  is  known  that  the  deflection  is  zero  at  the  supports, 
given  by  the  coordinates: 
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The  evaluation  of  Eq.  (210)  for  these  five  conditions 
gives  five  equations  and  five  unknowns  (Sp, 

A,  B,  C)  which  are  then  solved. 

The  same  procedure  is  followed  for  thi  deter¬ 
mination  of  and  S2p  and,  thus,  the  values  of 

the  Youngfs  moduli,  E  and  E  ,  and  the  Poisson-type 

x  y 

constants,  M  and  M  ,  are  known.  Thus,  the  ortho- 
xy  yx 

tropic  stiffnesses,  from  Eqs.  (20)  and  (21), 
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Another  scheme  for  determining  the  bending  rigidity 
is  presented  by  Stroud  (Ref.  335)  and  slightly 
modified  by  this  writer  in  an  unpublished  report. 

The  rectangular  plate  is  loaded  such  that  the  region 
between  supports  is  in  pure  bending  (Figure  45  )• 
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Figure  45 

Test  for  Bending  Stiffness  (II) 


For  the  shape  of  the  deflected  curve, 
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Differentiating , 
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Thus,  the  bending  moment,  from  Eq.  (17), 
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The  radius  of  curvature  between  points  i  and  j  is 
found  from  the  geometry: 
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Substituting  Eq.  (213)  into  Eq.  (212),  writing  the 
moment  in  terms  of  the  load,  and  solving  for  the 


stiffness 
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The  same  procedure  is  followed  for  D^.. 


8.2  Twisting  Stiffness 

The  loading  condition  required  for  the  deter¬ 
mination  of  the  twist  term,  D^,  is  that  of  uniformly 
distributed  edge  twisting  couples.  This  is  accomplished 
rather  simply  through  the  set-up  shown  in  Figure  46, 
which  is  described  by  several  writers,  among  them 
Timoshenko  and  Woinowsky-Krieger  (Ref.  97),  Hoppman 
et  al  (Refs.  329-332,  340),  Adams  (Ref.  29),  Stroud 
(Ref.  335),  and  Tsai  and  Springer  (Ref.  33&). 

As  is  illustrated  in  Figure  46,  a  square  plate 
is  subjected  to  concentrated  loading  at  one  pair  of 
diagonally  opposed  corners  and  supported  by  point 
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supports  (perhaps  ball  bearings)  at  the  other  two 
corners. 


£  PoimT  soppoir 
O  POINT  tOAP 


Figure  46 

Test  for  Twisting  Stiffness 


The  theory  given  by  Timoshenko  and  Woinowsky-Krieger 
(Ref.  97)  shows  that  the  concentrated  loads,  P/2, 
acting  at  each  of  the  unsupported  corners,  produce  a 
torsional  moment,  M  ,  per  unit  length  along  each  side 
of 
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The  deformation  is  given  by 
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Since  the  supports  are  at  the  same  height  and  the 
deflection  at  the  supports  is  zero  (w  =  0) ,  the  con¬ 
stants  A,  B,  C  can  be  determined.  Only  one  measure¬ 
ment  ,  at 


x  =  0,  y  =  0 

is  required  to  solve  Eq.  (216)  for  s66‘ 
the  relationships, 

,  i 
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are  used  to  find  the  twisting  stiffness. 

The  same  testing  procedure  can  be  used  for  a  more 

direct  method  of  calculation.  Since  the  edges  of  the 

twisted  plate  are  straight,  the  plate  is  in  a  condition 

of  pure  twist  with  no  shears  Q  and  Q  present.  Thus, 

x  y 

if  the  deflection,  5,  of  one  of  the  unsupported  corners 
(a/2,  -a/2)  or  (-a/2 ,  a/2)  is  measured,  its  relationship 
to  the  twist  is  given  by  the  simple  relationship, 

.  (.117) 

or 
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From  Eq.  (19),  the  torsional  rigidity, 
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From  this  and  Eqs.  (217)  and  (215), 
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This  is  substituted  into 


7  Co  * 


to  get  the  torsional  stiffness. 

Great  care  must  be  exercised  in  performing  the 
twist  tests.  It  must  be  remembered  that  the  edges  are 
to  be  kept  straight.  In  testing  patterns  such  as 
corrugations  or  dimples,  care  must  be  taken  to  allow 
each  cell  of  the  edge  cross-section  in  the  xz-  and 
yz-  planes  to  rotate  the  same  amount.  The  satisfaction 
of  these  two  requirements  may  require  stiffening  of  the 
edges  in  such  a  way  that  the  torsional  rigidity  of  the 
stiffeners  is  very  low  while  their  bending  rigidity  is 
high  (thin  plates),  and  perhaps  filling  the  edge  cells 


in  such  a  way  that  they  can  be  attached  to  the 
stiffeners. 

Becker  and  Gerard  (Ref.  323)  suggest  a  somewhat 
different  experimental  technique  in  that  a  canti¬ 
levered  plate  is  twisted.  But,  of  course,  the  twisted 
plate  has  the  same  shape  as  that  discussed  above  and 
the  method  of  calculating  the  constant  is  the  same. 

8.3  Shear  Stiffness 

When  the  shear  deformations  of  the  orthotropic 

plate  are  accounted  for  (Paragraph  1.8),  the  shear 

stiffnesses,  D  and  D  ,  must  be  determined.  Libove 

,  x  y 

and  Batdorf  (Ref.  155)  and  Libove  and  Hubka  (Ref.  156) 
discuss  a  method  for  this  calculation. 

The  transverse  shfsar  stiffness,  D  ,  is  found  by 
loading  the  plate  as  a  beam  by  a  uniformly  distributed 
load  (Fig.  47). 
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Figure  47 

Test  for  Shear  Stiffness 


The  bending  moment  at  any  point, 


The  rate  of  change  of  transverse  shear, 


(219) 
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The  curvature,  3w/3x*  ,  is  determined  from  deflection 
measurements  along  the  x-axis.  Then,  the  shear  stiff¬ 
ness  is  found  from: 


29. 


3x 

av  +  m, 

ax'-  15, 


y 


b 

^  J_  (P£*  .  px*\ 

*t>b  V  Z  2  / 


— r  S' 

in  the  same  manner. 


is  obviously  determined 


(221) 
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CONCLUSION 


The  expressions  given  in  this  report  should  enable 
the  analyst  to  calculate  the  stresses  and  deformations 
of  flat  plates  of  various  configurations  --  sandwich, 
corrugated,  rib-reinforced,  integrally  stiffened  --  due 
to  static  loading  conditions.  Various  loading  possi¬ 
bilities,  in  combination  with  a  variety  of  end  support 
conditions,  are  described. 

Perhaps  the  greater  contribution  of  this  report  is 
the  inclusion  of  the  most  complete  list  of  references 
of  work  in  the  area  of  orthotropic  plates.  The  bibli¬ 
ography  of  Part  9  is  divided  according  to  subject  matter 
in  such  a  manner  that  a  researcher  in  this  field  will 
find  it  convenient  in  locating  publications  pertinent 
to  his  work. 

A  second  report,  which  will  be  more  suitable  for 
design,  will  be  published  in  the  near  future.  In  this 
future  publication,  the  author  will  present  the  infor¬ 
mation  concerning  orthotropic  plates  in  such  a  manner 
that  the  designer  will  be  able  to  make  rational  com¬ 
parisons  between  the  various  configurations  of  stiffened 
plates,  based  upon  some  standard  criterion,  such  as 
weight  or  cost.  This  work  is  presently  taking  place, 
and  it  will  be  presented  in  a  report  titled  "Orthogonal.1 
Stiffened  Plates:  A  Design  Handbook." 
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